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During  the  course  of  the  "year  the>followlng  studies  ,•  v*w* 
c  juduutwJtr  'J)ody  wave  studied  of  special  effects  related  to  path 
of  propagation  and  station  site,  the  focal  mechanism  of  earth¬ 
quakes,  model  studies  cf  the  relation  between  cavity  size  and 
decoupling  of  explosive  sources,  the  determination  of  depth 
of  focus,  and  the  phase  velocity  of  surface  waves. ^Some  of 
these  studies  have  been  completed  and  are  published  or  sub¬ 
mitted  for  publication.  Others  are  still  under  investigation 
or  are  as  yet  in  the  process  of  preparation  for  publication. 
Accordingly,  this  summary  is  divided  into  two  sections  corre¬ 
sponding  to  the  completion  status  of  the  work. 

A.  Work  in  progress  or  nearing  completion. 

1.  Travel  Times  and  Amplitudes  of  T,ong  Period  Body 
Shear  Waves  from  Nuclear  Explosions 

The  nuclear  events  Half-Beak  and  Greeley,  which  took  place 
at  the  Nevada  Test  Site  on  30  June  and  20  December  1966, respec¬ 
tively,  generated  relatively  large  amplitude  long-period  shear 
waves.  In  particular  they  provided  data  for  determining  trav¬ 
el  times  and  amplitudes  of  long-period  body  shear  (S)  waves; 
until  this  time  no  nuclear  explosions  had  produced  useful  data 
of  this  type.  To  the  best  of  our  knowledge  the  travel  times 
and  amplitudes  included  in  this  report  are  the  first  of  this 
kind  for  explosion-generated  S  waves. 
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Seismograms  which  supplied  the  data  discussed  in  this  re¬ 
port  were  obtained  from  WWSS,  LRSM,  and  Saint  Louis  University 
network  stations.^  Additional  seismograms  from  the  LRSM  and 
Canadian  network  stations  have  been  requested  by  us;  they  will 
be  analyzed  in  ths  coming  months  as  a  continuation  of  the  work 
contained  in  this  report. 

For  the  purpose  of  determining  the  velocity  structure  of 
the  earth's  mantle  the  use  of  S  waves  has  one  decided  advan¬ 
tage  over  the  use  of  P  waves.  At  those  distances  where  the 
travel-time  curve  is  multi-valued,  the  arrivals  corresponding 
to  the  different  branches  can  often  be  recognized  by  a  change 
in  polarization  of  the  S  motion.  This  distinctive  feature  does 
not  apply  to  P  waves,  for  they  are  all  linearly  polarized.  The 
identification  of  later  ?  arrivals  can  be  accomplished  only  by 
detecting  changes  in  amplitude  or  period,  which  are  often 
difficult  or  impossible  to  discern.  Thus  the  later  arrivals 
of  the  multi-branched  travel-time  curves  are  more  readily  de¬ 
termined  for  the  S  waves  than  for  their  P  wave  counterparts. 

Identification  of  changes  in  polarization  and  of  onsets 
of  S  waves  is  best  accomplished  by  employing  particle  motion 
diagrams,  as  was  demonstrated  for  earthquake-generated  S  waves 
by  Ibrahim  and  Nuttli  (1967).  This  method  of  analysis  does 
present  practical  problems,  because  large  amounts  of  time  are 
required  for  digitizing  the  seismograms  and  constructing  the 
particle  motion  diagrams. 

Travel  Times.  The  travel-time  data  are  presented  in  the 
form  of  a  reduced  time-distance  graph  (Fig.  1).  The  solid 
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line  is  a  computed  travel-time  curve,  based  on  the  S  velocity 
model  of  Ibrahim  and  Nuctli  (1967)*  for  a  surface  focus.  Cir¬ 
cles  are  used  to  represent  the  arrival  times  of  WWSS  or  Saint 
Louis  University  network  stations,  and  crosses  for  the  LRSM 
stations . 

An  immediate  conclusion  which  follows  from  an  inspection 
of  Pig.  1  is  that  all  the  data  points  do  not  lie  on  the  curves. 
Much  of  the  scatter  is  reduced  if  one  applies  a  station  anomaly 
correction  as  determined  by  Doyle  and  Hales  (1967)  for  stations 
in  the  United  States.  For  example,  between  19°  and  22°  the  cir¬ 
cles  correspond  to  s  arrivals  at  ROL,  DBQ,  FLO,  SLM,  and  OXF; 
these  are  stations  in  the  central  United  States  for  which  Doyle 
and  Hales  found  the  S  arrivals  to  be  4  or  5  seconds  too  early. 
After  correction  the  points  would  fall  on  branches  4,  3, and  6, 
the  branches  of  the  largest  amplitude  arrivals  in  this  distance 
interval . 

Generally  it  is  not  possible  to  identify  arrivals  corre¬ 
sponding  to  all  the  branches  of  the  travel-time  curve  at  a  par¬ 
ticular  station.  There  are  a  variety  of  reasons  for  this,  but 
the  principal  ones  are  l)  that  at  certain  distances  the  ampli¬ 
tude  of  a  wave  of  a  particular  branch  can  become  very  small, 
and  2)  that  the  low  frequency  of  the  S  wave  motion  usually 
does  not  permit  a  time  resolution  of  events  whose  beginnings 
are  separated  by  less  than  ten  seconds,  especially  if  the  ear¬ 
lier  of  the  two  arrivals  has  a  larger  amplitude  than  the  later 
one. 

Table  1  contains  the  time  difference  between  the  first 
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TABLE  1.  TRAVEL  TIME  DIFFERENCES 
(Adopted  Times  Minus  Jeffreys-Bullen  Times) 
For  First  S  Arrivals  and  A  Surface  Focus 


Branch 

Time 

Difference 

Branch 

Time 

Difference 

5° 

1 

-  4.5  sec. 

28° 

6 

+5.1  sec 

6 

1 

-  5.6 

29 

6 

+  4.4 

7 

1 

-  6.6 

30 

6 

+  3»t* 

8 

1 

-  7.6 

31 

6 

+  3.0 

9 

1 

-  8.4 

32 

6 

+  2.5 

10 

1 

-  9.2 

33 

6 

+  2.2 

n 

1 

-  9.7 

34 

6 

+  2.0 

12 

1 

-10.0 

35 

6 

+  1.8 

13 

1 

-10.4 

36 

6 

+  1.2 

14 

1 

-10.3 

37 

6 

+  1.0 

15 

1 

-  9.9 

38 

6 

+  0.8 

16 

1 

-  9.8 

39 

6 

+  0.8 

17 

1 

-  8.9 

40 

6 

+  0.7 

18 

1 

-  8.1 

4i 

6 

+  0.7 

19 

1 

-  6.8 

42 

6 

+  0,8 

20 

1 

-  5.2 

43 

6  - 

+  0.9 

21 

1 

-  1.8 

44 

6 

+  1.3 

22 

4 

-  1.7 

45 

6 

+  1.6 

23 

4 

-  0.5 

46 

6 

+  1.9 

24 

4 

+  1.3 

47 

6 

+  2.0 

25 

4 

+  3.4 

48 

6 

+  2.5 

26 

4 

+  5.3 

49 

6 

+  2.8 

27 

6 

+  5.8 

50 

6 

+  3.4 

6 


arrival  of  our  travel-time  curve  and  the  Jeffreys-Bullen 
tables.  The  time  differences  are  negative  out  to  23°,  and 
positive  beyond  that.  Out  to  15°  the  large  negative  time 
differences  pr  bably  can  be  explained  by  the  fact  that  the 
S  wave  (branch  1,  usually  called  Sn)  has  a  small  amplitude 
compared  to  the  higher-mode  surface  waves  which  begin  impul¬ 
sively  five  to  ten  seconds  later  at  these  distances.  This 
higher-mode  surface  wave  most  likely  is  mistakenly  Identified 
as  S  in  the  majority  of  cases.  Prom  15°  to  25°  the  Jeffreys- 
Bullen  tables  likely  refer  to  branch  2  rather  than  branch  1, 
which  is  the  earlier  but  which  is  characterised  by  very  Bmall 
amplitudes . 

Amplitudes.  Amplitude  and  period  measurements  were  made 
of  the  first  .half-cycle  of  the  S  motion,  and  where  poBBible  of 
the  next  half-cycle  as  well.  Pig.  2  demonstrates  the  method 
used  to  obtain  these  quantities.  The  zero  to  peak  amplitude 
is  A],;  the  corresponding  period  is  4tp  where  t-^  is  the  time 
from  the  beginning  of  the  phase  to  the  first  peak.  Th/,  ampli¬ 
tude  of  the  second  half-cycle  of  motion  1b  taken  to  be  the 
average  of  the  amplitudes  of  the  first  and  second  peaks, namely 

(Ai  +  A2)/2J  the  corresponding  period  is  2(t2  -  t^),  where 
t2  -  t^  is  the  time  difference  between  the  peaks  of  motion. 

Tables  2  and  3  contain  the  amplitude  data  for  Half-Beak 
and  Greeley,  respectively.  Several  important  observations  or 
conclusions  follow  from  the  data  included  in  these  tables. 
First,  both  explosions  generated  a  large  amount  of  SH  energy, 
comparable  in  size  to  the  SV  energy.  Second,  the  amplitude 
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TABLE  ?.  AMPLITUDES  AT  S  WAVES  FOR  HALF-BEAK. 


Amplitude 

Branch  Station  Distance  Component  (zero-peak) 


Amplitude 

Period  (peak-peak)  Period 
2 


1 


2 


3 


4 


TUC 

6,764° 

SH 

0.1 

microns  6  seconds 

0.1 

microns  6  seconds 

R 

0.2 

6 

0.2 

6 

LUB 

12.342 

SH 

•0.3 

6 

Z 

0.2 

6 

RGSD 

12.419 

SH 

0.3 

22 

- 

- 

MHT 

15.605 

SH 

0.2 

6 

0.6 

8 

PGBC 

17.229 

Z 

0.1 

17 

RKON 

21.033 

SH 

0.4 

20 

• 

•m 

R 

0.4 

20 

- 

- 

PGBC 

17.229 

SH 

0.2 

23 

0.7 

20 

R 

0.1 

23 

_ 

Z 

0.1 

23 

0.1 

20 

SIBC 

18.983 

SH 

0.4 

22 

0 .6 

16 

R 

0.2 

22 

0.3 

16 

ROL 

19.353 

R 

0.2 

8 

0.4 

6 

DBQ 

20.287 

SH 

0.4 

16 

• 

R 

0.7 

8 

• 

Z 

0.7 

8 

• 

mm 

SLM 

20.569 

SH 

0.9 

10 

1.0 

10 

RKON 

21.033 

Z 

0.1 

24 

_ 

OXF 

21.930 

SH 

0.8 

7 

R 

1.2 

6 

Z 

0.8 

5 

SHA 

24.205 

SH 

0.8 

8 

R 

- 

- 

1.7 

1.0 

CRNB 

15.484 

Z 

0.3 

22 

0.2 

22 

MHT 

15.605 

SH 

0.7 

8 

0.6 

8 

Z 

0.6 

8 

0.7 

8 

JPAT 

15.631 

SH 

0.2 

16 

0.5 

16 

R 

0.1 

16 

0.1 

16 

Z 

0.1 

16 

0.1 

16 

KCMO 

17.089 

SH 

0.2 

22 

0.2 

12 

Z 

0.1 

22 

0.1 

12 

SIBC 

18.963 

Z 

0.4 

40 

ROL 

19.353 

SH 

0.4 

8 

Z 

0.3 

8 

0.6 

6 

SLM 

20.569 

Z 

0.5 

4 

0.4 

4 

BLO 

23.427 

R 

1.0 

28 

2.9 

12 

AX2AL 

25.055 

R 

0.1 

16 

Z 

0.1 

16 

BLA 

28.450 

SH 

0.2 

24 

0.2 

24 

R 

1.4 

28 

LPS 

33.274 

R 

0.5 

12 

0.4 

12 

COL 

33.287 

R 

0.3 

10 
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Table  d.  (Cont'd) 


Amplitude 

Branch  Station  Distance  Component  (zero-peak)  Period 


Amplituae 
(peak-peak)  Period 
2 


5 


6 


BLO 

23.427 

SH 

0.4  microns 

SHA 

24.205 

SH 

1.2 

HNME 

26.629 

Z 

0.1 

ROL 

19.353 

SH 

0.1 

SLM 

20.569 

R 

2.0 

Z 

0.6 

OXP 

21.930 

SH 

2.1 

R 

1.7 

BLO 

23.427 

SH 

0.3 

SHA 

24.205 

SH 

0.5 

AX2AL 

25.055 

R 

0.6 

BLA 

28.450 

SH 

0.4 

R 

0.5 

BE  PL 

29.751 

SH 

0.2 

R 

0.1 

LPS 

33.274 

R 

0.3 

COL 

33.287 

SH 

0.1 

R 

0.1 

36.629 

Z 

0.) 

HNME 

Z 

0.x 

8  seconds  0.7  microns  6  seconds 


17 

- 

mm 

32 

mm 

4 

10 

2.0 

8 

8 

1.0 

8 

15 

1.9 

17 

15 

0.8 

17 

12 

- 

- 

11 

- 

- 

16 

0.6 

16 

6 

0.4 

6 

6 

- 

- 

20 

0.3 

12 

20 

0.2 

12 

8 

0.4 

11 

10 

- 

- 

10 

- 

- 

10 

- 

- 

22 

- 

- 
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TABLE  3.  AMPLITUDES  OP  S  WAVES  FOR  GREELEY 


Amplitude 

Amplitude 

Branch 

Station 

Distance  Component 

(zero-peak) 

Period 

(peak-peak)  Period 

2 

1 

LON 

10.257° 
15 .368 

SH 

_ 

0,2  microns 

1  9  seconds 

JCT 

Z 

- 

- 

0.1 

12 

DAL 

16.664 

z 

0.4  microns 

6  sec. 

0.4 

10 

2 

DAL 

16.664 

R 

0.7 

4 

BLA 

28.539 

SH 

- 

• 

0.7 

9 

3 

DAL 

16.664 

SH 

1.7 

4 

mm 

FLO 

20.542 

SH 

3  A 

12 

3.0 

12 

SLM 

20.658 

SH 

5.9 

12 

4.6 

10  ' 

OXF 

22.018 

SH 

3.5 

24 

4.5 

22 

R 

2.3 

24 

2.0 

22 

Z 

1.1 

24 

0  9 

22 

4 

JCT 

15.368 

Z 

— 

0.2 

12 

SLM 

20.658 

R 

0.6 

8 

0.8 

8 

OXF 

22.018 

R 

0.4 

9 

0.5 

9 

AAM 

25.300 

SH 

0.7 

7 

0.4 

10 

BLA 

28.539 

R 

«a» 

- 

1.0 

11 

CMC 

30.591 

SH 

0.7 

5 

0.6 

8 

COL 

33.264 

SH 

0.4 

10 

- 

5 

SHA 

24.290 

SH 

2.2 

14 

2.4 

14 

R 

4.6 

14 

5.5 

14 

AAM 

25.500 

SH 

0.6 

mm 

- 

0.8 

11 

R 

7 

0.9 

11 

Z 

0.2 

7 

0.4 

11 

BLA 

28.539 

SH 

— 

1.2 

14 

R 

2.4 

16 

1.8 

14 

Z 

0.8 

16 

0.9 

14 

CMC 

30.591 

SH 

0.9 

16 

- 

- 

Z 

0.7 

16 

- 

— 

6 

FLO 

20.542 

R 

1.2 

7 

1.0 

10 

SLM 

20.658 

R 

3.8 

8 

SHA 

A  r 

?4.290 

SH 

1.6 

a 

0.9 

;  -1 

.1 

, » »•» 

bLA 

.  .  *  ■ 1  » 

28.539 

SH 

O 

11 

0.9 

8 

CMC 

30*591 

SH 

0.9 

10 

0.8 

8 

* 

GEO 

30.840 

SH 

0.3 

10 

- 

- 

COL 

33.264 

SH 

0.4 

11 

- 

- 

R 

0.2 

11 

- 

mm 

Z 

0.3 

11 

- 

- 

GDH 

46.103 

SH 

0.2 

8 

- 

- 

R 

0.3 

8 

- 

- 

Z 

0.3 

8 

- 

- 

SJG 

47.706 

SH 

0.2 

8 

0.2 

8 

R 

0.1 

8 

0.1 

8 
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of  the  waves  of  branch  1  from  5°  to  20°  Is  about  one  order  of 

magnitude  less  than  those  of  branches  3,  4,  5 ,  and  6  from 
o  o 

about  20  to  30  .  The  amplitude  of  t».-  waves  of  branch  2  is 
also  comparatively  small,  but  not  as  small  as  those  of  branch 
1.  Third,  the  waves  of  branch  3  appear  to  have  the  largest 
amplitude.  This  seems  surprising  in  terms  of  our  velocity 
model,  for  they  correspond  to  a  narrow  beam  of  rays  leaving 
the  focus  at  angles  of  incidence  between  4l°  and  45°j  these 
rays  bottom  in  the  mantle  at  a  depth  of  about  380  km.  Our 
calculati  ..2  of  the  transfer  function  for  a  crustal  layer 
show  that  for  these  particular  frequencies  and  angles  of  in¬ 
cidence  a  large  amplification  effect  of  the  crust  on  the  sur¬ 
face  motion  may  occur  for  the  SV  component  of  motion.  How¬ 
ever,  this  would  not  explain  the  large  SH  motion. 

One  point  should  be  kept  in  mind  when  one  attempts  to 
analyze  these  amplitude  data.  Amplitudes  and  periods  are  dif¬ 
ficult  to  measure  at  those  distances  at  which  the  travel-time 
curve  is  multi-valued.  One  may  be  dealing  with  the  superposi¬ 
tion  of  two  or  more  waves  arriving  at  approximately  the  same 
time,  rather  than  a  single  wave. 
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2.  The  Use  of  Amplitude  Data  In 
Focal  Mechanism  Studies 

Commonly  focal  mechanism  studies  are  based  upon  the  signs 
of  the  first  motion  of  the  P  waves.  Because  the  distribution 
of  selsmographlc  stations  seldom  Is  adequate  to  determine 
uniquely  the  orientation  of  the  nodal  planes  (or  of  the  x  and 
y  axes,  the  lines  normal  to  these  planes),  Stauder  and  his 
students  have  Included  S  wave  polarization  data  In  their  focal 
mechanism  studies.  In  this  manner  they  have  greatly  Improved 
the  resolution  of  the  nodal  plant  a .  However,  In  an  appreciable 
number  of  earthquakes  It  Is  still  possible  to  assign  a  range 
of  values  to  one  of  the  parameters  of  the  focal  mechanism  solu¬ 
tion,  e.g.  the  strike  of  the  y  axis.  The  purpose  of  this  study 
is  to  determine  If  amplitude  data  can  provide  an  additional 
constraint  on  the  focal  mechanism  solution,  so  that  a  more  near¬ 
ly  unique  solution  can  be  established. 

Theory.  Our  model  of  the  focal  mechanism  is  a  douole 
couple  with  moment,  the  so-called  Type  II-  force  model.  The 
displacements  at  large  distances  from  the  source  in  an  infinite, 
homogeneous  medium  are  (Stauuer,  1962): 


4  IT  a^ 


K'  (t  -  R/a) 


(1) 


c 

4  IT  b3R3 


K'  (t  -  R/b) 


(2) 


where  x,y,z  are  coordinates  referred  to  the  system  of  forces, 
R  is  the  distance  from  a  point  (x,y,z)  to  the  focus,  a  and  b 
are  velocities  of  the  P  and  S  waves,  respectively,  t  is  time, 
and  K(t )  is  a  function  of  the  amplitude  and  time  history  at 

e 


12 


the  source.  The  x  and  y  axes  have  the  direction  of  the  forceo 
of  the  couple  acting  at  the  origin,  and  z  is  orthogonal  to 
them  so  as  to  form  a  right-handed  system. 

In  addition  we  require  a  geographic  coordinate  system, 
x,  y,  z  j  also  with  the  focus  as  origin.  The  coordinates  are 
related  by  the  equations 

x  -  OSc  *  +  £>xy  +  Yxz 
y  “  «y  *  *  (Sr  y  +  yy  * 

where  _ 

x  *  cos  AZes  sin  lh 

y  ■  sin  AZea  sin  lh 
z  -  cos  ih  * 

In  these  equations  ih  is  the  angle  of  incidence  at  the  focus 
and  AZea  is  the  azimuth  from  the  epicenter  to  the  point  (x,  y,  z). 

The  direction  cosines  are  given  by 
v*  -  cos  AZ„  cos  PL 

X  A  A 

f>x  =  sin  AZX  cos  PLX 

-  sin  PLX 

oty  »  cos  AZy  cos  PLy 

^y  =  sin  AZy  cos  PLy 

Yy  “  BlR  PLy  - 

where  AZ^  and  PLX  are  the  azxmuth  and  plunge  of  the  x-axis,  re¬ 
spectively.  Similarly  AZy  and  PLy  are  the  azimuth  and  plunge 
of  the  y-axis.  The  condition  that  the  x  and  y  axes  be  perpen- 
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dlcular  leads  to  the  equation 

tan  PLX  tan  PLy  =  -  cos  (AZ*  -  AZy). 

Finally,  c,  which  appears  in  eq.  (2),  is  given  by 

(  0^  y  +  vXy  x)  y  -  (£xy+£vx)x 
sin  ih 

The  assumptions  made  in  deriving  eqs.  (l)  and  (2)  do  not 
apply  to  the  real  earth.  We  must  take  account  of  the  fact  that 
it  is  a  finite,  bounded  non- homogeneous  sphere  rather  than  an 
infinite,  homogeneous  space.  This  produces  an  additional  de¬ 
pendence  of  amplitude  on  epicentral  distance,  owing  both  to 
geometrical  spreading  and  to  absorption.  The  effect  of  geomet¬ 
rical  spreading,  which  can  be  calculated,  is  larger  than  that 
of  absorption  for  the  low-frequency  body  waves  with  which  we 
are  concerned. 

The  free  surface  and  the  crustal  layers  cause  the  motion 
at  a  point  on  the  surface  of  the  earth  to  differ  from  that  at 
a  neighboring  point  just  below  the  ..rust.  This  effect,  which 
depends  both  on  the  angle  of  incidence  and  the .frequency  content 
of  the  incident  wave  motion,  is  much  more  severe  for  the  SV 
component  of  motion  than  for  P  and  SH.  Accordingly  we  restrict 
our  amplitude  studies  to  P  and  SH. 

For  a  specified  epicentral  distance  all  the  terms  in  eqs. 
(l)  and  (2)  except  for  x, y,  and  c  will  be  constant.  Thus,  if 
the  observed  P  wave  amplitudes  are  divided  by  the  product  xy, 
a  plot  of  these  adjusted  amplitudes  versus  epicentral  distance 
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should  show  a  more-or-less  smooth,  continuous  ohange  In  ampli¬ 
tude  with  distance,  as  for  a  uniform  radiating  source.  A  large 
amount  of  scatter  of  the  data,  particularly  for  stations  near 
a  nodal  line.  Indicates  that  the  focal  mechanism  parameters 
should  be  revised.  In  this  manner  P  wave  amplituder  can  be 
used  to  refine  fooal  mechanism  solutions.  For  SH  waves  one 
adjusts  th«  observed  amplitude  by  dividing  It  by  the  value  of 
c. 


Examples .  We  consider  here  two  earthquakes  whose  focal 
mechanism  solutions  have  been  given  by  Stauder  and  Bollinger 
(1966,  1965)*  The  first  Is  the  Kodiak  Island  earthquake  of 
February  6,  1964,  a  foreshock  of  the  great  Alaskan  earthquake 
of  the  following  month.  The  seoond  Is  an  earthquake  which 
occurred  In  .the  Kurile  Islands  on  June  28,  1963.  Table  1  lists 
the  hypocentral  coordinates  and  the  focal  mechanism  solutions 
for  these  earthquakes. 

Table  1 


Date 

Location 

Latitude 

Longitude 

Origin  Time 

Focal  Depth 

Magnitude 


Feb.  6,  3964 
Kodiak  Islands 
55.TN 
155. 8W 
13-07-25 
33  km 
7 


June  28,  1963 
Kurile  Islands 
46. 5N 


153.2E 
21-55-39 
33  km 
6-3/4 


Stauder  &  Stauder  St 

Bollinger  Revised  Bollinger  Revised 


AZX 

336° 

339° 

162° 

1540 

PLX 

9 

12 

74 

76 

AZy 

130.5 

118 

296 

320 

pv 

80 

75.5 

12 

13.: 
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For  the  Kodiak  Island  earthquake  we  read  the  amplitudes 
of  the  P  waves  from  the  long-period  seismograms  of  fifty  VWSS 
stations.  Total  ground  amplitudes  were  obtained  by  taking  the 
square  root  of  the  sum  of  the  squares  of  the  vertical  and  hori¬ 
zontal  components.  In  some  cases  where  horizontal  data  were 
lacking  the  total  amplitude  was  obtained  by  dividing  the  ver¬ 
tical  component  by  the  cosine  of  the  apparent  angle  of  inci¬ 
dence.  The  amplitudes,  uncorrected  for  focal  mechanism,  are 
plotted  as  a  function  of  distance  in  Fig.  1.  The  data  scatter 
widely  in  this  figure,  and  show  no  observable  dependence  on 
eplcentral  distance . 

Fig.  2  illustrates  the  marked  reduction  in  the  scatter  of 
amplitudes  which  occurs  after  the  amplitudes  have  been  adjusted 
according  to.  the  focal  mechanism  solution  of  Stauder  and  Bol¬ 
linger.  However,  three  points  have  large  amplitudes  and  plot 
off-scale.  In  addition,  six  points  (indicated  by  broken  x's 
in  the  figure)  have  the  wrong  sense  of  motion.  The  size  of 
the  symbols  is  an  indication  of  the  distance  of  a  station  from 
a  nodal  line;  small  symbols  are  for  stations  near  nodal  lines, 
and  large  ones  for  stations  more  removed  from  them.  The  ad¬ 
just*,  i  amplitudes  of  stations  near  nodal  lines  are  most  sen¬ 
sitive  to  the  proper  focal  mechanism  solution. 

By  a  Judicious  trial  and  error  procedure  we  have  modi¬ 
fied  the  values  of  the  azimuth  and  plunge  of  the  x  and  y  axes 
to  further  reduce  the  scatter  in  the  amplitude  data.  Our  best 
solution,  from  a  set  of  eighty  considered  solutions,  is  given 
as  the  revised  solution  in  Table  1.  Adjusted  amplitudes  using 
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this  solution  are  presented  in  Pig.  3.  For  this  solution  no 
stations  have  an  Incorrect  sense  of  first  P  motion,  and  the 
adjusted  amplitudes  for  stations  near  nodal  lines  are  reason¬ 
ably  similar  to  those  for  stations  far  removed  from  them  but 
at  the  same  epicentral  distance.  The  solid  line  in  this  fig¬ 
ure  is  a  theoretical  curve  showing  the  expected  decrease  in 
amplitude  owing  to  geometrical  spreading,  calculated  from  the 
Jeffreys-Bullen  travel  time  tables.  Because  absolute  ampli¬ 
tudes  for  this  theoretical  curve  are  unknown,  it  was  arbltrar- 

o 

lly  made  to  fit  our  data  at  80  .  Its  remarkably  close  agree¬ 
ment  with  the  observed  data  leads  u*s  to  conclude  that  geo¬ 
metrical  spreading  is  the  principal  cause  of  amplitude  atten¬ 
uation  for  long-period  P  waves. 

Fig.  4  presents  the  SH  amplitudes  adjusted  for  the  focal 
mechanism,  based  upon  the  revised  solution  given  in  Table  1. 
The  scatter  is  reasonably  small,  less  than  for  any  other  fo¬ 
cal  mechanism  solution  which  we  considered.  We  speculate  that 
the  increase  in  amplitude  which  is  observed  beyond  70°  is 
caused  by  the  superposition  of  ScS  motion  with  that  of  S. 

Fig.  3  contains  the  adjusted  P  wave  amplitude  data  for 
the  Kurile  Islands  earthquake.  Although  there  is  more  scatter 
of  these  data  than  for  the  Kodiak  Island  earthquake,  the 
scatter  for  this  solution  is  significantly  less  than  for  most 
of  the  eighty  other  solutions  which  we  considered.  The  same 
statement  applies  for  the  adjusted  SH  amplitudes,  which  are 
given  in  Fig.  6. 
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Future  Work.  We  intend  to  continue  our  studies  of  the 
type  reported  here,  paying  particular  attention  to  earthquakes 
located  in  regions  where  a  knowledge  of  the  focal  mechanism 
parameters  car  assist  in  a  better  understanding  of  the  geology 
of  these  regions. 

From  our  present  studies  it  is  ppparsnt  that  no  meaning¬ 
ful  amplitude-distance  studies  can  be  undertaken  without  cor¬ 
recting  for  the  earthquake  radiation  pattern,  and  that  even 
small  errors  in  the  focal  mechanism  parameters  can  produce 
large  errors  and  large  scatter  in  the  adjusted  amplitude  data. 
We  now  are  lr  a  position  to  obtain  much  better  amplitude  atten¬ 
uation  data,  which  was  the  primary  purpose  for  our  undertaking 
these  amplitude  studies. 
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3.  Azimuthal  Variation  of  the  P  Wave  Particle  Motion 

from  the  Great  Circle  Path. 

P  wave  amplitudes  contain  information  about  the  s*~ucture 
of  the  crust  and  upper  mantle  at  the  recording  site.  This  in¬ 
formation  may  be  used  to  supplement  that  obtained  from  refrac¬ 
tion  and  surface  wave  dispersion  studies  or,  when  these  are 
lacking,  to  assign  preliminary  values  to  the  crust  and  upper 
mantle  parameters.  In  this  study  we  are  employing  the  depar¬ 
ture  of  the  direction  of  the  horizontal  component  of  the  P 
wave  motion  from  the  great  circle  path,  and  In  particular  Its 
dependence  on  the  azimuth  of  the  approaching  wave,  to  Infer 
the  earth  structure  at  selected  WWSS  stations. 

In  cases  where  the  distance  from  the  epicenter  Is  con¬ 
siderable,  a  longitudinal  wave  may  be  regarded  as  plane  where 
It  emerges.  When  a  longitudinal  wave  strikes  the  earth's 
crust,  if  such  boundaries  are  horizontal,  the  refracted  or 
reflected  waves  which  are  formed  will  always  be  polarized  in 
such  a  way  that  the  displacements  will  remain  In  the  vertical 
station  to  epicenter  plane.  It  follows  that  It  is  possible 
to  determine  the  azimuth  towards  the  epicenter .from  the  sta¬ 
tion,  using  the  horizontal  components  of  the  vibrations  In 
tne  P  wave.  That  Is,  the  azimuth  &  is  given  by 

Y  V 

tan  _ H 

*n  ve 

where  Y0  and  are  amplitudes  measured  on  the  E-W  and  N-S 
seismograms  during  the  first  Impulse  of  the  longitudinal  wave, 
and  Ve  ami  Vn  are  the  dynamic  magnifications  of  the  seismo¬ 
graphs.  An  Incorrect  solution  for  &  may  result  from  errors 
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In  the  determination  of  the  instrumental  constants  cr  In  case 
the  Inner  refractive  boundaries  are  not  horizontal.  The  latter 
constitutes  a  matter  of  Interest  In  the  present  investigation 
of  the  peculiarities  of  the  structure  of  the  earth. 

The  simplest  earth  model  which  can  produce  an  azimuth 
anomaly  is  a  single  dipping  discontinuity  surface.  For  such 
a  case  the  azimuth  anomaly  curve  is  psuedo-sinusoidal,  with 
zero  anomaly  in  the  dip  direction.  Fig.  1  presents  theoreti¬ 
cal  curves  for  a  plane  dipping  at  15°,  with  a  velocity  of  5.6 
km/sec  above  and  7.75  km/sec  below.  The  solid  line  is  for  an 
angle  of  incidence  of  25°  (epicentral  distance  of  about  50°) 
and  the  dashed  line  for  angles  of  incidence  of  45°  to  60° 
(epicentral  distance  of  about  20°  to  25°). 

Experimental  studies.  The  experimental  work  consists  in 
measuring  the  half  amplitudes  and  the  peak-to-peak  amplitudes, 
as  well  as  the  periods,  of  the  first  motion  on  the  N-S  and  E-W 
component  long-period  seismograms,  using  a  digitizing  machine. 
By  applying  the  known  constants  of  the  instruments  these  mea¬ 
sured  values  are  converted  into  ground  motion,  to  the  nearest 
one-tenth  of  a  micron.  These  values  are  used  to  determine  an 
azimuth,  which  is  called  the  observed  azimuth.  A  calculated 
azimuth  is  obtained  using  a  formula  of  spherical  trigonometry, 
where  the  coordinates  of  the  epicenter  and  seismograph  station 
are  known.  The  azimuth  anomaly  is  defined  as  the  observed 
value  minus  the  calculated  value. 

Following  is  a  brief  discussion  of  the  features  of  the 
azimuth  anomaly  curves  constructed  for  certain  selected  WWSS 


stations • 
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ALQ  (Figs.  2  and  3)*  For  Albuquerque  It  was  possible 

to  divide  the  observations  Into  three  groups ,  depending  on 

the  signs  of  the  anomaly.  In  Fig.  2  It  can  be  seen  that  with 

the  exception  of  two  points  all  of  the  anomalies  between  0° 

and  180°  are  negative.  Positive  anomalies  are  found  between 

180°  and  300°,  and  negative  beyond  300°.  The  maximum  value 

o 

of  the  anomaly  Is  of  the  order  of  3  on  the  positive  side  and 
o 

10  on  the  negative  side.  Fig.  2  Is  based  on  the  first  half 
amplitude,  and  Fig.  3  on  the  peak-to-peak  amplitude  data. 

Both  show  the  same  trends,  although  there  are  fewer  data  for 
the  latter  figure  because  In  many  cases  the  P  motion  Is  In¬ 
terrupted  by  other  arrivals  before  one  full  cycle  Is  completed. 

BKS  (Fig.  4)  In  the  case  of  the  Berkeley  data  the  obser¬ 
vations  are  more  scattered,  which  lb  probably  Indicative  of 
the  more  complex  crustal  and  upper  mantle  structure  at  that 

station.  Most  of  the  points  between  0°  and  160°  are  eharae- 

o  o 

terlzed  by  a  negative-  anomaly,  while  those  from  160  to  316 

o 

are  positive  and  those  beyond  316  are  negative.  The  maximum 

o 

value  of  the  anomaly  Is  about  8  on  both  the  positive  and 
negative  sides . 

MDS  (Fig.  3)  We  have  fewer  data  for  Madison  because  the 

seismograms  are  often  disturbed  by  long-period  noise.  We  have 

only  two  observations  at  azimuths  0°  to  110°,  for  which  both 

anomalies  are  small  and  positive.  With  one  exception  all  the 

o  o 

anomalies  between  azimuths  of  130  and  330  are  negative,  the 
amplitude  being  as  large  as  24°.  These  large  negative  anom¬ 
alies  are  surprising  If  one  imagines  the  crustal  structure  to 
be  simple  and  composed  of  nearly  horizontal  layers. 


AZIMUTH  ANOMALY  W  DEGREE  AZIMUTH  ANOMALY  IN  DEGREE 
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OXF  (Fig.  6)  It  was  possible  to  obtain  many  seismograms 
showing  impulsive  P  waves  at  Oxford.  The  anomalies  are  nega¬ 
tive  between  0°  and  40°,  positive  between  4ov‘  and  200°,  nega¬ 
tive  between  200°  and  300°,  and  positive  beyond  300°.  The  maxi¬ 
mum  amplitude  of  the  anomaly  is  of  the  order  of  12°  on  the 

o 

positive  side  and  6  on  the  negative  side. 

It  is  interesting  to  compare  Fig.  6  with  Fig.  2  for  Al¬ 
buquerque.  One  is  almost  the  mirror  image  of  the  other.  That 
is ,  at  azimuths  where  ALQ  has  positive  anomalies,  OXF  has  nega¬ 
tive,  and  vice  versa.  The  shapes  and  amplitudes  are  also  almost 
the  same.  This  indicates  that  the  structural  features  at  the 
two  stations  are  eimilar  and  have  the  same  strike  direction, 
but  that  the  dip  directions  differ  by  180°. 

TUC  (Fig.  7)  A  good  azimuthal  distribution  of  the  data 
could  not  be  obtained  in  the  case  of  Tucson.  Most  of  the  ob¬ 
servations  are  between  the  azimuths  of  120°  to  180°  and  300° 

o  o 

to  330  .  Most  of  the  anomalies  are  positive  between  12C  and 

o  o 

236  and  negative  beyond.  Maximum  amplitudes  are  20  on  the 

positive  side  and  15°  on  the  negative  side. 

Future  Work.  It  is  obvious  that  the  theoretical  curves 
for  a  single  dlpplr.g  layer  do  not  satisfy  our  observations. 

We  shall  consider  other  earth  models,  and  derive  formulas  for 
the  azimuth  anomaly  produced  by  them  in  an  attempt  to  get  a 
better  match  to  our  data. 

We  shall  apply  the  same  techniques  to  the  long-period 
LASA  data  from  Montana.  Because  of  the  unusually  low  seis¬ 
micity  on  a  global  scale  for  this  year  we  have  only  found  two 
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e-rthquakes  in  the  first  seven  months  of  1967  which  provide 
satisfactory  LASA  data.  As  we  wait  for  sufficient  LASA  data 
to  accumulate  we  shall  study  other  WWSS  stations  in  the  U.  S. 
for  the  periods  1964-1966. 

4.  P  Time  Delays  and  Continental  Uplift. 

Introduction.  The  velocity  versus  depth  profile  under¬ 
lying  the  travel  time  tables  of  seismic  waves  by  Jeffreys  and 
Bullen  (1948)  proves  to  be  of  remarkable  consistence  all  over 
the  world.  In  the  past,  deviations  of  the  travel  times, thought 
not  without  importance,  barely  exceeded  the  observational  accu¬ 
racy  of  most  selsmographlc  stations.  This  contributes  certainly 
to  the  fact  that  the  Jeffreys-Bulien  travel-time  tables  became 
the  most  used  ones,  at  least  as  far  as  teleseisms  are  concerned. 

Since  two  parameters,  the  eploentral  distance  and  the 
focal  depth  appear  to  be  of  predominant  importance  for  the 
travel  time,  the  structure  of  the  earth  is  mainly  concentric. 

Perturbations  of  the  travel  times  are  of  regional  importance. 
Their  causes  are  located  in  the  crustal  and  upper  mantle  struc¬ 
ture  underneath  the  selsmographlc  station  or  the  source.  The 
investigation  of  the  perturbations  was  facilitated  in  the  last 
years  by  the  increased  quality  of  the  instrumentation  of  seis- 
mographlc  observatories.  The  present  study  was  possible  because 
the  travel  times  at  some  selsmographlc  stations  are  measured  and 
reported  with  an  accuracy  of  0.1  sec. 


The  objective  of  the  study  is  to  Investigate  the  possible 
correlation  between  the  continental  uplift  and  the  P-tlme  delays 
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Fennoscandia  was  chosen  because: 

1.  I*  ie  covered  by  a  relatively  dense  network  of  Beis- 
mographic  observatories, 

2.  It  has  a  seismographic  service  of  high  standard,  and 

3*  The  evidence  of  continental  uplift  is  well  estab¬ 
lished  by  geologic,  geodetic,  and  mareographio  obser¬ 
vations  made  over  a  period  of  decades  and  even  cen¬ 
turies  . 

The  method  applied  here  was  stimulated  by  a  recent  inves¬ 
tigation  by  Bolt  and  Nuttli  (1966).  These  authors  found  that 
the  P-time  delays  at  twelve  seismographic  stations  in  central 
and  northern  California,  if  related  wO  a  reference  station, 
Berkeley,  vary  cyclically  with  the  azimuth  from  the  station  to 
the  epicenter.  In  the  present  study  the  P-time  delays  at  15 
Scandinavian  stations  were  investigated.  Uppsala  and  UraeS  were 
used  as  reference  stations.  This  produced  28  sets  of  relative 
P-time  delays. 

Observational  Material.  In  the  time  from  January  1,  1963 
to  December  31,  1965*  15  seismographic  stations  were  operational 
in  Norway,  Sweden,  and  Finland,  for  which  station  bulletins  were 
available.  The  length  of  the  interval  guaranteed  a  sufficiently 
large  sample  of  data,  "he  beginning  of  the  interval  w*s  chosen 
to  coincide  with  the  time  when  meat  of  the  stations  started  to 
report;  the  arrival  times  of  P-waves  with  an  accuracy  of  0.1  sec. 
The  names,  coordinates,  elevations  above  sea  Itvel  of  the  sta¬ 
tions,  and  intervals  for  which  the  bulletins  were  available, 
are  given  in  Table  1.  The  location  of  the  stations  is  seen  in 
Pig.  1. 
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The  arrival  times  of  P  waves  from  more  than  300  earth¬ 
quakes  have  been  extracted  from  the  bulletins  of  the  15  sta¬ 
tions.  For  the  most  part  the  earthquakes  were  in  the  epi- 
central  distance  range  between  30°  and  100°,  and  were  selec¬ 
ted  bo  as  to  cover  possibly  all  azimuths.  Obviously,  the 
distribution  of  seismic  regions  limits  an  even  coverage  of 
the  azimuths.  For  the  same  reason,  the  epicentral  distances 
of  the  earthquakes  studied  change  with  azimuths,  and  the 
assumption  had  to  be  made  that  P-time  delays  are  independent 
of  the  epicentral  distance. 

The  P-time  delays  found  at  each  station  over  a  period  of 
3  years  may  be  expected  to  scatter  about  the  calculated  val¬ 
ues,  especially  since  the  readings  were  performed  in  course 
of  routine  work.  The  relative  P-time  delays,  involving  read¬ 
ings  at  two  stations,  will  scatter  even  more.  However  it  is 
clear  that  with  a  sufficiently  large  amount  of  readings  sig¬ 
nificant  conclusions  can  be  drawn  with  regard  to  the  relative 
P-time  delays. 

We  found  the  greatest  uniformity  of  readings  at  the 
Swedish  stations  Uppsala  (UPP),  Klruna  (KIR),  Skalstugan 
(SKA),  Qoteborg  (QOT),  Uraefi  (UME),  Karlskrona  (KLS) .  Two 
of  these  station^  Uppsala  and  Ume8,  were  taken  as  reference 
stations  for  all  the  remaining  stations.  Further  reasons 
for  choosing  them  as  reference  stations  are: 

1.  Neglecting  the  array  station  Llll&hammer  (LIL), 
Uppsala  and  UmeS  recorded  more  earthquakes  and 
provided  more  travel  times  than  any  other  station. 
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2.  Uppsala  and  especially  Umea  are  situated  at  places 
of  most  Intensive  continental  uplift,  granting 
eventually  a  uniformity  of  the  structure  beneath 
the  stations,  and  at  the  same  time  close  to  the 
center  of  the  Scandinavian  stations  involved  (see 
Pig.  1). 
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Lillehammer  Is  the  only  array  station  here.  It  was 
operated  during,  and  bulletins  were  Issued  for,  28  out  of 
36  months  Investigated  (see  Table  1).  The  station  produced, 
however,  only  147  usable  time  differences,  as  compared  with 
the  maximum  of  297  at  Uppsala  and  UmeS. 

Relative  P-tlme  Delays.  The  following  quantity  was 
computed  for  all  earthquakes  recorded  at  any  of  the  14  sta¬ 
tion  (ST)  and  the  reference  (REP)  station  (Uppsala  or  UmeS): 

RPD  -  (Tobfl  -  Tcal)gT  -  (TQba  -  (i) 

where  Tobs  -  the  observed  travel  time,  and 

Tcft^  -  the  travel  time  calculated  from  Jeffreys- 
Bullen  Tables. 

For  both  Tobfl  and  TQal  ,  the  origin  time,  epicentral  coordin¬ 
ates,  and  focal  depth  of  the  event  are  accepted  as  given  by 
the  USCOS  In  their  Seismologlcal  Bulletin. 

A  relative  P-time  delay  so  constructed  represents  a 
travel  time  anomaly.  Its  cause  Is  present  In  that  portion 
of  the  two  ray  paths,  where  the  latter  diverge  most,  i.e.  in 
the  upper  mantle  and  crust,  below  the  two  stations  Involved. 

Whenever  the  RPD  exceeded  +  3  sec.,  it  was  rejected  as 
unreliable  and  containing  u  systematic  .rror.  The  total 
number  of  time  differences  produced  in  that  way  for  each 
station  is  given  in  Table  1. 
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The  RPD  were  assigned  a  weight  of  0.5  if: 

1.  The  arrival  time  was  reported  with  an  accuracy  of 
only  1.0  sec. 

2.  The  RPD  was  diverging  by  **1  sec.  from  the  delays 
at  neighboring  azimuths. 

3.  The  distance  from  t"ie  epicenter  to  a  station  or  a 
reference  station  became  smaller  than  30°. 

If  two  of  the  conditions  applied  to  an  RPD,  the  weight  assigned 
was  0.25;  If  all  three  applied  -  the  figure  was  dropped.  All 
other  RPD's  were  given  the  weight  1.0  .  The  number  of  RPD's 
with  the  weights  1.0  ,  0.5  ,  0.25  for  the  stations  are  given 
in  Table  1. 

The  RPD  were  plotted  as  function  of  an  angle  9  equal 
to  the  arithmetic  average  of  the  azimuths  frc^i  the  station 
and  reference  station  to  the  epicenter.  As  Illustration, 

Pig.  2  gives  the  RPD's  as  function  of  the  azimuth,  for  o£te- 
borg  versus  Uppsala. 
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The  RPD  for  data  secs  were  approximated  In  the  least 
square  sense  by  the  equation 

RPD  -  A  +  B  sin(©  +  € )  (2) 

where  A,  B,  £  are  constants  characteristic  of  each  station. 
They  were  calculated  for  14  stations,  if  referred  to  Upp¬ 
sala,  and  again,  if  referred  to  UmeS.  Table  2  contains  the 
corresponding  figures.  Most  remarkable  is  that  with  Uppsala 
as  reference  station  all  14  A-values  are  positive  with  an 
average  of  0.47;  with  UmeS  as  reference  station  13  stations 
attained  a  positive  A-value,  the  average  amounting  to  O.53. 

This  indicates  a  certain  similarity  of  the  structures  below 
Uppsala  and  UmeS.  It  might  be  noted  here  that  the  A-values 
in  all  approximations  are  closest  to  zero  for  Ume8,  Uppsala, 
and  Karlskrona  (see  Fig.  l),  all  three  stations  delineating 
a  roughly  NS  trending  axis. 

The  RPD  were  correlated  with  a  sinusoidal  curve  shifted 
along  the  axis  of  azimuths  at  intervals  of  2°.  FOr  the  data 
shown  in  Fig.  2,  the  correlation  coefficient  amounts  to  0.52 
for  a  Dag  £  being  equal  to  that  found  from  the  least  square 
approximation  within  less  than  10  degrees.  The  correlation 
with  272  data  pairs  is  then  significant  at  the  70  per  cent 
confidence  level. 

Orientation  of  the  RPD  Approximations.  The  least  square 
approximation  of  the  RPD's  by  a  sinusoidal  curve  has  an  appre¬ 
ciable  standard  deviation.  In  the  example  shown  on  Fig.  2, 
the  standard  deviation  is  0.82.  However,  as  mentioned  earlier. 
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Table  2 


Coefficients  A,  B,  and  £  in  equation 
RPD  *  A  +  B  sin(©  +£  ) 

Reference  stations: 


Uppsala  (UPP) 

Um  e£ 

(UME) 

A 

B 

£ 

A 

B 

£ 

UPP 

-0.04 

0.52 

95.61 

KIR 

0.38 

0.49 

-106.51 

0.31 

0.18 

157.65 

SKA 

0.55 

0.44 

-134.99 

0.45 

0.46 

140.65 

GOT 

0.45 

0.53 

-171.41 

0.50 

0.79 

148.49 

UME 

0.04 

0.52 

-  84.39 

KLS 

0.13 

0.48 

-146.33 

0.03 

0.54 

140.02 

HEL 

1.02 

0.23 

173.18 

0.84 

0.57 

113.65 

NUR 

0.45 

0.24 

-  58.45 

0.51 

0.25 

99.60 

SOD 

0.41 

0.24 

-  10.83 

0.38 

0.48 

61.38 

KJN 

0.45 

1.02 

-  45.10 

0.56 

0.64 

-  28.85 

KEV 

0.34 

0.23 

-22.48 

0.41 

0.23 

85.03 

KON 

0.75 

0.33 

127.09 

0.97 

0.8l 

132.76 

TRO 

0.47 

0.24 

-174.49 

0.71 

0.60 

162.43 

KRK 

0.55 

0.37 

-  19.24 

0.75 

0.10 

90.51 

LIL 

0.53 

0.46 

-175.68 

0.55 

0.73 

137.89 
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the  RPD  were  obtained  from  routine  measurements  made  for  two 
stations  over  a  time  interval  of  3  years,  the  stations  having 
Independent  recording  systems  and  clocks.  A  thorough  analysis 
of  variance  renders, for  every  particular  station  taken  in¬ 
dividually,  only  a  low  reliability  of  the  approximation.  Nev¬ 
ertheless  for  each  of  the  reference  stations,  Uppsala  and 
UmeS,  14  approximations  are  available,  and  we  may  look  for  a 
systematic  behavior  of  all  the  approximations  even  though  no 
statistical  model  is  available  which  would  permit  the  estima¬ 
tion  of  the  confidence  level  of  conclusions  drawn  from  the 
Interrelated  approximations. 

On  Figure  3  the  arrows  represent  the  azimuths  of  the  max¬ 
imum  positive  RPD's  for  all  14  stations,  referred  to  Uppsala. 
Figure  4  shoyrs  the  corresponding  directions  with  Umefi  as  ref¬ 
erence  station.  The  arrows  are  assigned  to  the  midpoints  be¬ 
tween  the  particular  station  and  the  reference  station.  It 
can  be  seen  immediately  that  the  arrows  have  a  tendency  to  be 
directed  away  from  the  region  of  most  intensive  continental 
uplift,  in  the  direction  perpendicular  to  the  lines  of  equal 
rate  of  uplift  (compare  Figs.  3  and  4  with  Fig.  l). 

Of  particular  interest  is  the  bend  of  the  directions  at 
points  situated  to  the  W,  N,  and  E  of  UmeS,  around  the  region 
of  most  intensive  uplift  so  as  to  point  away  from  it. 

This  fact,  together  with  the  positive  A-values  pointed 
out  earlier,  indicates  that  seismic  velocities  below  the  ref¬ 
erence  stations  are  higher  than  below  the  other  stations, 
situated  at  places  of  less  intensive  continental  uplift,  as 
will  be  discussed  with  greater  detail  below. 


FIG.  4 
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Without  discussing  now  the  exceptions  as  seen  particu¬ 
larly  for  Helsinki  (HEL)  on  Fig.  3,  and  for  Helsinki  and 
Nurmi jarvi  (NUR)  on  Fig.  4,  we  will  propose  an  interpreta¬ 
tion  of  the  general  features  of  the  RPD  as  function  of  azi¬ 
muth. 

Travel  Times  as  Function  of  Azimuth.  Let  us  assume  two 
layers  in  contact  with  the  outer  boundaries  being  plane, 
horizontal  and  parallel  to  each  other,  and  the  plane  inter¬ 
face  to  dip  at  an  angle  £  (see  Fig.  5A).  The  velocities  of 
P-waves  are  vA  and  v0  in  the  lower  and  upper  layer  respective¬ 
ly.  The  total  thickness  of  the  layers  is  H. 

A  P-wave  arriving  from  an  azimuth  measured  from  the 
dip  direction  will  be  refracted  at  the  lower  boundary  and  at 
the  Interface,  suffering  in  the  latter  case  a  deviation  <p  from 
the  original  azimuth.  The  deviation  can  be  expressed  as 

sin 'I' 

m  arc  tan 

with  i1  *  angle  of  refraction  at  the  lower  bou^lary;  this  angle 
can  be  found  if  the  angle  of  incidence  at  the  lower  boundary 
and  the  velocities  on  both  sides  of  it  are  known.  It  will  be 
constant  for  a  particular  eplcentral  distance  (see  Fig.  5B). 

The  travel  time  across  the  two  layers  will  be  obviously 
a  cyclical  function  of  the  azimuth.  It  can  be  expressed  by 
the  following  formula: 
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where  h1  is  tjie  depth  of  the  refraction  point  at  the  dipping 
interface  btaow  the  surface,  ig  and  i^  are  the  angles  of  in¬ 
cidence  and  refraction  at  the  interface  respectively,  and  ^3 
expressed  byt 


arc  sin 


sin  i. 


— - -  i  i...  sraii.  i-.gmaaKgaa— — — ^ 

cjs  i^  .  y  tan2!^  +  2  tan  i^.tan  S  cos  vp  +  tan* 


The  RPD,  as  calculated  from  the  observations,  are  the  differ¬ 
ences  of  two  travel  tiroes  of  the  form  (2).  They  will  attain 

»  o  _  o 

extremal  values  f or  -  0  ,  or  180  .  Assuming  the  two  stations 

and  the  epicenter  to  be  situated  along  the  dip  direction  of  the 
interface  at  a  relative  distance  of  500  km,  and  assuming  vA  , 
vfi  ,  and  i^  to  be  V600  km/sec,  4500  km/sec,  1°,  and  18.5°  re¬ 
spectively,  the  extremal  RPD  c_.n  be  estimated  to  be  of  the 
order  of  -0.75  and  0.75  sec.  The  figures  are  compatible  with 
the  amplitude  of  the  RPD's  observed  in  our  investigation.  The 
calculated  range  of  amplitudes  depends  on  a  number  of  param¬ 
eters  which  unfortunately  cannot  be  determined  uniquely  from 
the  c" nervations. 

Two  remarks  concerning  the  formula  (2)  are  appropriate 
here.  The  thickness  H  considered  does  not  coincide  neces¬ 
sarily  with  any  marked  boundary,  e.g.  the  Mohorovicic  discon¬ 
tinuity.  It  defines  only  a  level  above  which  a  travel  time 
anomaly  exists  along  the  ray  paths  to  the  two  stations.  The 
dip  of  the  interface  does  not  necessarily  mean  a  real, 
dipping  boundary.  It  may  be  understood  as  a  measure  of  a 
horizontal  velocity  Inhomogeneity.  This  means  that  the  medium 
may  have  a  velocity  gradient  in  the  dip  direction  of  the  inter- 
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face,  different  from  zero.  This  model  seems  most  relevant  to 
the  RPD's  as  observed  in  Scandinavia. 

In  Fig.  5C  the  station  (ST),  reference  station  (R),  and 
the  epicenter  are  assumed  to  lie  in  a  vertical  plane  coin¬ 
ciding  with  the  dip  direction  of  the  interface.  The  reference 
station  (R)  is  situated  "-lose  to  the  outcrop  of  the  Interface, 
so  that  the  rays  arriving  at  (R)  travel  always  with  the  veloc¬ 
ity  vA  assumed  to  be  larger  than  vfi. 


In  this  case  the  RPD's  will  be  always  positive.  The 
maximum  value  is  obtained  for  waves  traveling  in  the  updlp. 
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the  minimum  value  -  for  waves  traveling  in  the  downdip  direc¬ 
tion  changing  cyclically  for  the  remaining  azimuths. 

It  can  be  seen  that  essentially  the  same  will  apply  if  a 
horizontal  velocity  gradient  exists  in  the  layer,  and  if  the 
reference  station  1j  situated  in  the  direction  of  the  hori¬ 
zontal  velocity  gradient  with  respect  to  the  station.  And 
opposite,  if  the  reference  station  is  situated  away  f:om  the 
station,  in  the  direction  opposite  to  the  horizontal  velocity 
gradient,  the  average  of  the  RPD's  for  all  azimuths  will  be 
negative. 

We  shall  apply  this  model  to  the  RPD's  as  observed  ir, 
Scandinavia.  The  A- values  in  formula  (l)  were  found  to  be 
positive  for  Uppsala  and  Umei  as  reference  stations.  This 
implies  that,  the  reference  stations  a: '  situated  in  the  direc¬ 
tion  of  a  horizontal  velocity  gradient.  Since  the  reference 
stations  are  surrounded  by  14  stations  each,  we  conolude  that 
underneath  the  reference  stations  the  velocities  of  P  waves 
are  highest. 

Interpretation.  As  was  mentioned  earlier,  the  refer¬ 
ence  stations  Uppsala  and  UmeS  are  situated  at  places  of  most 
intensive  continental  uplift.  Let  us  assume  here  the  hypoth¬ 
esis  that  continental  uplift  in  Fennoscandia  is  an  evidence 
of  lsostatlc  readjustment  of  the  crust.  If  lsostasy  alone 
would  be  involved,  we  would  expect  the  horizontal  velocity 
gradient  in  the  crust  and  upper  mantle  to  be  directed  away 
from  the  center  of  most  intensive  uplift,  since  the  latter 
coincides  with  the  region  where  the  lighter,  crustal  material 
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was  mostly  submerged  during  the  glacier  period. 

Our  Investigation  Indicates  the  velocity  gradient  to  be 
directed  towards  the  center  of  uplift.  This  does  not  contra¬ 
dict  the  lsostasy  as  mechanism  of  the  uplift,  but  since  isos- 
tasy  cannot  explain  the  finding,  an  additional  mechanism  needs 
to  be  Introduced. 

The  most  obvious  one  seems  to  be  the  compression  of  the 
layers  In  consequence  of  loading  with  the  Icecap,  the  compres¬ 
sion  having  been  most  Intensive  In  the  region  of  most  Inten¬ 
sive  present  uplift.  During  glaciation  the  material  In  the 
crust  and  upper  mantle  can  be  supposed  to  have  been  compressed 
elastically  In  a  time-independent  (perfect  elasticity)  and  a 
time-dependent  way,  the  latter  having  a  recoverable  and/or  an 
Irrecoverable  component.  Unloading  caused  a  dilatation  In  the 
same  order,  the  time-dependent  dilatation  eventually  still  go¬ 
ing  on.  Compression  Is  known  to  cause  In  general  an  Increase 
of  seismic  velocities,  and  the  region  mostly  suppressed,  re¬ 
vealing  now  the  most  intensive  uplift,  should  then  have  higher 
seismic  velocities  if  compared  with  the  surrounding. 

Conclusions. 

1.  A  cyclical  variation  of  relative  P-time  delays  at  14 

selsmographlc  stations  in  Scandinavia  Is  found,  with 
o 

Uppsala  and  Umea  as  reference  stations. 

2.  Twenty-eight  sets  of  relative  P-time  delays  are  con¬ 
structed  as  function  of  azimuth  from  the  station  to 
epicenter,  and  the  observational  data  are  approx¬ 
imated  by  the  first  two  terms  of  a  Fourier  series 
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expansion.  The  constant  term  is  found  to  be  positive, 
and  the  maximum  positive  delay  is  found  for  azimuths 
pointing  away  from  the  region  of  maximum  continental 
uplift . 

3.  The  behavior  of  relative  P-tlme  delays  Is  Interpreted 
In  terms  of  a  horizontal  velocity  Inhomogeneity  In 
the  crust  and  upper  mantle,  such  that  higher  veloci¬ 
ties  are  found  in  the  region  of  maximum  continental 
uplift . 

4.  The  velocity  Inhomogeneity  does  not  exclude  lsostasy 
as  mechanism  of  the  continental  uplift;  an  additional 
mechanism  of  the  uplift,  however,  must  be  Introduced 
to  explain  the  findings. 

5.  As  additional  mechanism  we  propose  the  following.  The 
strata  were  compressed  during  glaciation.  In  conse¬ 
quence  of  the  removal  of  he  icecap  a  dilatation  took 
place.  Both  compression  and  dilatation  had  a  time- 
independent  and  a  time-dependent  component.  Part  of 
the  compression  may  even  be  Irrecoverable.  In  regions 
of  most  Intensive  present  uplift  the  compressions  may 
be  supposed  to  have  been  heaviest,  and  even  now  higher 
sell  mlc  velocities  may  be  expected  there,  as  found 
Indeed  from  the  relative  P-tlme  delays.  Whether  the 
time-dependent  dilatation  is  still  going  on  Is  beyond 
the  method  applied  In  the  present  investigation. 
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5.  The  Mechanism  of  the  Earthquakes  of  the  Rat  Island 
Earthquake  Sequence  of  4  February  1965. 

A  magnitude  7  3/4  (PAS,BRK),  8-8fc  (PAL),  7.5  (CQS-MC) 
earthquake  occurred  near  Amchitka  Island  on  4  February  1965, 
05-01-21.8  GMT.  The  USCOS  epicenter  has  coordinates  51.3°N, 
178. 6°E,  with  focal  depth  about  40  km  (restricted).  The 
major  shock  was  followed  by  a  strong  aftershock  sequence.  For 
example,  of  the  many  aftershocks  which  occurred  during  the 
first  24  days  following  the  main  shock,  750  were  large  enough 
to  be  located  by  the  routine  procedures  of  the  epicenter  de¬ 
termination  program. 

This  study  concerns  the  focal *mechanism  of  the  main  shock 
and  of  the  larger  aftershocks. 

Fig.  1  presents  a  map  showing  the  location  of  the  main 
shock,  the  extent  of  the  epicentral  region,  and  the  location 
of  the  larger  aftershocks  which  are  here  examined  in  detail. 
The  aftershock  epicenters  shown  on  the  map  are  those  of  the 
smaller  aftershocks  located  by  the  USCOS  for  the  first  three 
days  following  the  major  shock,  together  with  those  of  the 
larger  shocks  of  the  entire  sequence.  It  is  seen  that  the 
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major  shock  occurred  near  the  eastern  extremity  of  the  activ¬ 
ity,  that  the  aftershock  zone  extends  from  180°  to  170°E,  and 
that  the  activity  extends  over  a  zon?  about  200  km  wide  im¬ 
mediately  to  the  south  of  the  island  chain.  The  focal  depths 
given  by  the  USCGS  vary  from  25  to  40  km.  A  few  foci  lie  in 
the  range  of  focal  depth  40  to  50  km. 

The  earthquakes  here  examined  are  indicated  by  the  lar¬ 
ger  symbols  on  the  figure.  The  solid  symbols  represent  foci 
of  the  first  three  days.  Thereafter  there  were  no  aftershocks 
large  enough  to  yield  a  focal  mechanism  solution  until  30  March. 
The  open  symbols  represent  earthquakes  occurring  between 
30  March  1965  and  4  July  1966.  The  epicentral  coordinates, 
together  with  the  focal  mechanism  solutions,  for  all  these 
earthquakes  are  listed  in  Table  1. 

Focal  Mechanism  Solutions.  The  focal  mechanism  solutions 
have  been  determined  using  data  from  the  WWSSN  system.  In  the 
determinations  data  from  the  long -period  records  only  were  used. 
Both  direction  of  first  motion  of  P  and  the  polarization  of 
the  S  wave  were  used  in  determining  the  nodal  planes.  In  all 
cases  a  double  couple  mechanism  was  assumed;  no  evidence  favor¬ 
ing  or  requiring  any  other  type  mechanism  was  observed. 

Two  distinct  classes  of  foci  were  found.  An  example  of 
the  first  i3  indicated  in  Fig.  2a.  The  observable  P  wave  ra¬ 
diation  pattern  of  this  class  samples  only  two  first  motion 
quadrants  on  the  focal  sphere.  As  a  result,  only  one  nodal 
plane  of  the  F  wave  first  motion  is  determined.  On  the  assump¬ 
tion  of  the  double  couple  mechanism,  the  S  wave  polarization 
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Table  1.  Rat  Island  Earthquake  Sequence 

4  February  1965 


Plane  1 

Plane  2 

pate 

Time 

Lat. 

Lon*. 

h 

b 

(USCOS) 

Dip 
Direc . 

DIP 

Tip 

Dlrec 

.  Dl£ 

4  Feb  65 

05-01-28 

51. 3N 

178. 6E 

40 

6.0* 

150°E 

•-J 

00 

0 

- 

- 

n 

12-06-04 

52. 6N 

172. IE 

25 

5.8 

145°E 

83° 

29°E 

15° 

it 

14-18-28 

53.  ON 

171  OE 

30 

5.7 

l48°E 

83° 

0 

33  E 

17° 

5  Feb 

09-32-09 

52. 3N 

174. 3E 

.41 

5.9 

l4l°E 

75° 

325°B  15° 

11 

20-47-13 

51. 9N 

174. 6E 

35 

5.7 

0 

153  E 

77° 

310°E 

14° 

5  Feb 

04-02-53 

52.  IN 

175. 7E 

35 

5.9 

142°E 

75° 

359°E 

18° 

7  Feb 

02-17-09 

51. 4N 

173. 4E 

40 

6.0 

196°E 

53° 

58° 

45° 

it 

09-25-51 

51. 4n 

179. IE 

30 

5.3 

l4l°E 

74° 

34o°E 

16° 

30  Mar 

02-27-07 

50. 6N 

177. 9E 

51 

7 

194°E 

47° 

52°E 

50° 

23  May 

23-46-12 

52. 2N 

175. OE 

22 

6.1 

145°E 

75° 

12°e 

22° 

1  Oct 

08-52-04 

50.  IN 

178. 2E 

23 

6.3 

185°E 

50° 

40°E 

50° 

22  Nov 

20-25-13 

51. 4N 

179. 7W 

40 

5.9 

145°E 

70° 

354°E 

22° 

15  May  66 

14-46-06 

51. 5N 

178. 4W 

31 

5.8 

136°E 

61° 

349°E 

34° 

2  Jime 

03-27-53 

51.  IN 

176. OE 

41 

6.0 

0 

211  E 

51° 

50°E 

40° 

4  July 

18-33-36 

51. 7N 

179. 9E 

13 

6.2 

97°E 

85° 

187°E 

85° 

♦This  is  the  rob  value  given  by  the  USCOS  from  short-period 
body  waves.  Actual  magnitude  is  about  M  -  7.5.  Other 
shocks  here  listed  are  probably  in  the  range  6<M<7. 


data,  as  presented  on  the  right-hand  portion  of  the  figure, 
make  possible  a  complete  determination  of  both  planes. 

An  example  of  the  second  class  of  foci  is  presented  in 
Pig.  2b.  For  these  earthquakes  the  observable  P  wave  radia¬ 
tion  pattern  samples  only  a  single  quadrant.  Neither  nodal 
plane  is  determinable  on  the  basis  of  the  P  wave  data  alone, 
though  each  of  the  nodal  planes  is  restricted  to  a  range  of 
dip  in  the  neighborhood  of  45°,  +  15°.  In  this  case,  again 
on  the  assumption  of  a  double  couple  mechanism,  the  S  wave 
polarization  data  make  possible  a  determination  of  both  nodal 
planes.* 

The  orientation  of  the  nodal  planes  for  all  the  earth¬ 
quakes  studied  is  given  in  Table  1.  The  focal  mechanisms  are 
further  indicated,  together  with  their  spatial  relationship 
to  one  another,  in  Pig.  3.  The  shaded  portions  of  the  mechan¬ 
ism  diagrams  in  the  figure  correspond  to  rarefaction  first 
motion  quadrants. 

Discussion.  Several  comments  may  be  made  concerning  the 
focal  mechanism  of  these  Rat  Island  earthquakes. 

First,  there  is  a  distinct  spatial  classi¬ 
fication  of  the  foci.  One  group  consists  of  foci  located 

*  The  particular  configuration  in  question  is  one  which  corre¬ 
sponds  to  a  near  vertical  axis  of  greatest  compressive  stress 
fall  rarefaction  first  motions)  or  of  least  compressive  stress 
(all  compressions).  In  such  cases  a  rotation  in  the  strike 
of  the  nodal  plane  is  possible  with  no  great  difference  in  the 
agreement  with  the  S  wave  data.  Were  the  principal  axis  to  be 
perfectly  vertical,  the  solution  would  be  indeterminate,  save 
for  the  restriction  placed  on  the  dip  of  the  nodal  planes. 


near  the  inner  or  northern  edge  of  the  active  zone,  immed¬ 
iately  south  of  the  island  chain.  A  second  group  consists 
of  t!  e  foci  located  vnder  or  to  the  south  of  the  Aleutian 
Trench. 

The  first  group  Includes  the  initial  shock  and  ten  lar¬ 
ger  aftershocks.  The  focal  mechanisms  of  all  of  these  shocks, 
except  that  of  the  final  larger  shock,  that  of  4  July  1966, 
are  similar  to  that  shown  in  Fig.  2a.  This  means  that  each 
of  these  shocks  represents  a  continuation  of  episodic  fau.'it 
motion  of  one  and  the  same  character.  The  fault  plane  solu¬ 
tions  themselves  are  ambiguous  as  to  the  character  of  the 
motion:  either  of  the  nodal  planes  may  correspond  tc  the  fault 
plane.  If  the  steeply  dipping  plane  is  taken  as  the  fault, 
the  motion  is  reverse  faulting  along  a  series  of  faults  strik¬ 
ing  obliquely  to  the  trend  of  the  island  chain.  If  the  nearly 
horizontal  plane  is  taken  as  the  fault  plane,  the  motion  is 
thrust  faulting  in  character  along  a  plane  dipping  gently 
(about  10°  -  15°)  under  the  island  arc.  In  the  latter  hypoth¬ 
esis  the  island  or  landward  side  overthrusts  the  oceanic  side. 
As  indicated  by  the  arrows  in  Fig.  3,  the  thrust  block  moves 
toward  the  southeast.  The  motion  corresponds  to  a  pressure 

axis  approximately  normal  to  the  trend  of  the  arc  acting  up- 

o 

ward  from  the  ocean  side  at  an  angle  of  about  30  from  the 
horizontal . 

The  similarity  of  these  mechanisms  to  those  of  the  Alaska 
earthquake  of  23  March  1964  and  its  aftershock  sequence  (see 
Stauder  and  Bollinger,  1966)  is  noteworthy.  A  possible  special 
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significance  to  the  raeonani&m  of  focus  15,  that  of  4  July  1966, 
Is  noted  below. 

The  second  group  of  earthquakes,  those  with  foci  under 
the  trench,  are  all  similar  in  character  to  that  shown  in 
Pig.  2b.  The  pressure  axis  for  these  foci  is  oriented  nearly 
vertical,  whereas  the  axis  of  least  compressive  stress  is 
horizontal.  These  foci  apparently  do  not  differ  significantly 
in  focal  depth  from  those  of  the  first  group.  It  is  difficult 
to  correlate  the  stress  systems  of  the  two  groups. 

The  property  here  noted  for  the  earthquakes  with  foci  un¬ 
der  the  trench  (first  motions  all  rarefactions  at  teleselsmlc 
distances,  pressure  axis  vertical)  corresponds  to  a  property 
noted  previously  by  Evison  (1967)  for  earthquakes  under  the 
Aleutian  trench  as  opposed  to  those  with  foci  on  the  interior 
of  the  arc.  Foci  of  several  other  earthquakes  associated  with 
the  Aleutian  trench  which  we  have  studied  recently  have  this 
same  property. 

Secondly,  a  noteworthy  temporal  classification  of  the 
foci  suggests  itself.  Excluding  the  first  seve*.  hours  after 
the  main  shock,  during  which  time  the  traces  on  the  WWSSN 
records  were  too  tangled  to  interpret,  seven  aftershocks  large 
enough  for  a  focal  mechanism  solution  occurred.  There  followed 
a  lapse  of  fi.'ty  days  before  another  aftershock  of  comparable 
magnitude  occurred.  Thereafter,  seven  large  aftershocks  were 
placed  sporadically  over  the  next  fifteen  months. 

The  aftershocks  of  the  first  three  days  are  indicated  in 
Fig.  4a.  As  noted  above,  the  main  shock  occurred  near  the 
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eastern  extremity  of  the  eplcentral  zone.  The  larger  after¬ 
shocks  which  followed  were  distributed,  sequentially  in  time, 
at  the  opposite  extremity  of  the  active  zone  and  toward  the 
center.  The  next  to  last  aftershock  of  this  group  took  place 
under  the  trench,  and  the  l~«t  at  the  right  extremity  of  the 
zone.  The  mechanism  of  this  last  aftershock  differs  slightly, 
as  if  the  trend  of  the  faulting  curved  around  at  the  end  of 
the  fault. 

The  larger  aftershocks  of  the  later  group,  beginning 
30  March  1965,  concentrate  under  the  trench  and  to  the  east 
of  the  major  eplcentral  region,  along  the  continuation  of  r.he 
trend-  of  the  Aleutian  arc  in  the  Andreanof  Islands. 

The  last  shock  of  the  sequence,  that  of  4  July  1966,  is 
of  markedly  different  mechanism  from  the  others.  The  first 
motion  field  determines  well  all  four  first  motion  quadrants, 
and  the  faulting  motion  is  strike  slip  -  either  left  lateral 
on  a  north-south  fault,  or  right  lateral  on  an  east-west 
fault.  The  former  interpretation  is  here  preferred  as  more 
in  keeping  with  the  general  tectonic  movement  indicated  in 
Fig.  3.  Shortly  after  this  earthquake,  on  7  August  1966,  a 
large  earthquake  occurred  near  the  Fox  Islands  (50. 6N,  171.3°W) 
apparently  marking  the  initiation  of  a  new  phase  of  the  Aleu¬ 
tian  Island  activity. 

Conclusions.  Two  conclusions  of  significance  for  this 
Contract  follow  from  this  study.  First,  the  similarity  of 
mechanism  of  aftershocks  of  a  given  sequence  is  further 
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illustrated.  Second,  foci  radiating  uni-polar  first  motion 
radiation  fields  are  probably  restricted  to  certain  definite, 
and  relatively  narrowly  defined,  portions  of  a  tectonic  en¬ 
vironment  . 
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6.  Model  Studies  of  Decoupling  of  Explosions. 

During  the  last  period,  the  following  experiments  were 
performed . 

1.  Absolute  calibration  of  the  receiving  channels. 

2.  Decqupling  of  the  explosion. 

Receiving  channels  play  a  major  role  in  controlling  re¬ 
corded  wave  forms.  The  objective  of  this  calibration  is  to 
obtain  a  response  curve  in  absolute  units,  in  order  to  reduce 
the  displacement  of  the  wave  displayed  on  the  screen  of  the 
scope  to  the  displacement  at  the  surface  of  the  model.  The 
recording  channels  consist  of  a  capacitance  pick-up  and  a  60 
cycle  rejection  filter.  The  displacement  at  the  surface  of 
the  model  produces  a  change  in  capacity  of  the  capacitor, 
which  is  amplified,  filtered,  and  finally  displayed  on  the 
screen  of  the  scope.  In  order  to  calibrate  the  receiver  ab¬ 
solutely  a  source  of  mechanic?’’  motion  is  required.  A  piezo¬ 
electric  transducer  is  used  to  furnish  the  mechanical  motion. 
Transducers  were  sandwiched  between  two  circular  rods  of 
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Plexiglas  7/8"  in  diameter  and  equal  in  length.  Three  rods 
<*f  different  lengths  (10,  15>  and  20  cm)  were  used.  The  model 
is  suspended  by  a  fine  wire. 

In  this  experiment  a  sine  wave  frequency  oscillator  was 
employed.  The  signal  frequency  of  the  oscillator  vas  varied 
constantly  from  five  kilocycles  to  seventy-three  kilocycles 
per  second,  keeping  the  oscillator  output  amplitude  constant. 
The  output  voltage  of  the  oscillator  was  fifteen  volts  (RMS). 
It  was  found  that  the  amplitude  of  the  resulting  sinusoidal 
wave  pattern  on  the  oscilloscope  sc:  en  passed  through  a  ser¬ 
ies  of  maxima  and  minima,  and  the  amplitude  and  frequency 
values  of  these  maxima  were  measured.  The  maxima  occurre'  at 
resonant  frequencies  of  the  rods.  The  value  cf  these  reso¬ 
nant  frequencies  depends  on  the  elastic  properties  of  the 
rods  and  the  length  of  the  rods.  TH  resonant  fz'tquencies 
of  a  rod,  free  at  both  ends,  are  given  by 


in  which  f  «  frequency,  k  «  an  integer,  L  »  length  of  the 
rod,  and  C  »  phase  velocity  of  the  rod  material. 

In  the  configuration  in  question  the  transducer  is  al¬ 
ways  located  at  an  antinode  in  displacement;  therefore,  each 
rod  will  vibrate  as  if  free  at  both  ends.  For  the  even  num¬ 
ber  of  modes,  the  free  ends  of  the  rods  are  180°  out  of  phase 
with  the  input  of  the  transducer,  while  for  the  odd  number  of 
modes,  they  are  in  phase. 
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The  amplitude  recorded  at  the  free  end  of  the  rod  is 
the  incident  amplitude,  plus  the  multiple  reflections  from 
the  free  ends  of  the  rods.  The  possible  number  of  reflec¬ 
tions  is  infinite.  But  the  attenuation  in  plexiglas  is  high 
enough  to  diminish  the  amplitude  to  a  negligible  amount  after 
some  number  of  reflections.  To  correct  the  recorded  ampli¬ 
tude,  a  simplified  formula  is  derived  to  reduce  the  recorded 
amplitude  to  one  resulting  from  the  incident  wave  only. 

This  is  given  by 


Ap(R  +  pi*1 
R 


rK  e-2K«L 


in  which 

Xf  =  amplitude  at  the  ends  of  the  rod  plus  reflection 
A0+  -  output  amplitude  of  the  transducer 
R  *  reflection  coefficient  for  plexiglas  air  Interface 
L  «  length  of  the  rod 

*  absorption  coefficient  in  plexiglas 
k  -  irteger 

The  observed  amplitude  on  the  scope  ■  M5+  ■  As.  where  M  * 
magnification  factor 

"  “  (-1)K+1  RK  e-2K*L 


The  absorption  coefficient  is  obtained  from  another  experi¬ 
ment,  in  which  a  rectangular  pulse  generator  is  used  instead 
of  a  sine  wave  oscilator.  The  value  of  is  calculated  from 
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the  spectrum  of  the  pulses  reflected  at  the  free  ends  of  the 
rods. 

The  above  series  does  converge.  For-  calculation  of  mag¬ 
nification  factor  only  100  terms  of  the  series  were  included. 

The  amplitude  response  of  the  receiving  channels  obtained 
is  shown  in  Fig.  1. 

Phase  Response.  The  phase  response  of  the  receiving  chan¬ 
nels  was  obtained  as  follows.  The  scope  was  triggered  inter¬ 
nally  by  the  output  of  the  sine  wave  oscillator.  The  output 
of  the  scope  was  displayed  on  the  scope  on  one  of  the  beams. 

The  sine  wave  oscillator  also  applied  to  the  transducer  simul¬ 
taneously  and  the  output  of  the  channels  was  displayed  on  the 
upper  beam  of  the  scope.  It  is  assumed  that  there  is  no  de¬ 
lay  between  the  output  of  the  oscillator  and  the  output  of 
the  transducer.  The  phase  shift  of  the  sine  signal  propagated 
through  the  plexiglas  rod  and  recorded  by  the  receiving  unit, 
compared  with  the  output  of  the  oscillator  will  be  due  to  the 
receiving  channels. 

The  odd  modes  of  the  rods  are  180°  cut  of  phase  with  the 
even  mode.  To  measure  the  delay  introduced  by  the  receiving 
channels,  the  output  of  the  oscillator  is  compared  with  the 
output  */i*  the  receivers.  For  odd  mode  the  peaks  of  the  two 
outputs  arr  compared,  and  for  even  mode  the  peak  of  one  out¬ 
put  is  compared  with  the  trough  of  another  output.  Then  the 
delay  time  is  multiplied  by  the  angular  frequency  of  the  sine 
wave  to  obtain  the  phase  delay.  The  phase  response  curve  is 
shown  in  Fig.  2. 
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Decoupling .  This  experiment  was  performed  in  a  cubic 
foot  block  of  plaster  of  paris. 

The  decoupling  assembly  was  prepared  in  the  following 
manner.  Cylindrical  tubes  of  plexiglas  of  two  different 
diameters  50  mm  and  70  mm  were  used  to  catt  the  cylindrical 
cores.  Each  tube  consisted  of  two  parts ,  one  eight  centi- 
meters  in  length.  The  lower  part  of  the  specimen  was  casted 
first.  During  the  casting  of  the  lower  part  a  plexiglas 
sphere  was  placed  at  the  center  of  the  tube.  In  order  to 
prevent  the  sphere  from  sinking,  a  small  hole  about  one  milli¬ 
meter  in  diameter  was  made  through  the  sphere.  The  sphere  was 
held  at  the  center  by  a  thin  wire  wnich  ran  through  the  Bphero 
and  along  the  edge  of  the  tube.  After  putting  thin  aluminum 
foil  around  .the  sphere  and  over  the  lower  part,  the  upper 
part  of  the  core  was  cast.  The  advantage  of  this  procedure 
is  to  get  a  fine  contact  between  the  two  pieces  in  order  to 
cement  them  together  later  on.  After  five  minutes  of  setting, 
the  mold  is  removed,  the  two  ^leces  of  the  core  are  then  sep¬ 
arated  and  the  plexiglas  sphere  is  removed.  This  is  the  only 
way  that  could  produce  a  perfect  spherical  cavity. 

After  letting  the  core  specimen  dry  for  a  period  of  about 
two  weeks,  the  specimens  are  ready  for  use.  The  drying  period 
could  be  made  shorter  by  heating  the  sample  in  an  oven  for 
twenty-four  hours  with  50°  centigrade  temperature. 

The  charge  is  hung  inside  the  cavity  by  the  detonation 
wire  which  runs  through  the  charge.  A  fine  wire  was  used  in 
order  to  prevent  any  disturbance  of  the  pressure  inside  the 
cavity.  After  the  charge  was  assembled,  the  two  core  pieces 
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were  cemented  together  with  an  epoxy  cement.  The  samples 
were  left  overnight  to  let  the  cement  dry.  Weights  were  placed 
on  top  of  the  core  to  Join  the  two  pieces  firmly  together. 

The  core  can  then  be  fitted  Inside  the  block.  A  holj  Is 
drilled  with  power-driven  drilling  bits.  Two  sizes  of  holes 
were  drilled  fifty  millimeters  and  seventy  millimeters  in  di¬ 
ameter.  In  order  to  cement  the  core  to  the  hole,  a  soft  mix¬ 
ture  of  plaster  of  parls  and  water  was  poured  Inside  the  hole 
before  the  sample  was  lowered.  Immediately  after  pouring  the 
plaster  the  sample  was  pushed  Inside  the  hole.  The  core 
squeezed  the  plaster  out  along  the  gap  between  the  core  and 
the  hole  cavity,  and  the  plaster  filled  the  gap. 

The  advantages  of  this  technique  are  better  and  excel¬ 
lent  control,  of  the  charge  insi  e  the  cavity,  the  cavitieB 
are  perfect  In  shape  and  the  same  block  may  be  used  as  many 
times  as  one  wishes  without  spending  time  preparing  another 
block.  The  only  thing  one  does  1b  drill  the  hole  a^ain  and 
lower  another  specimen. 

To  establish  a  spherical  detonation  the  charge  should 
be  Initiated  at  the  center.  A  thin  tin  copper  wire  is  filed 
at  one  point  one  millimeter  In  length.  This  wire  Is  then 
passed  through  the  center  of  the  charge  with  the  filed  point 
at  the  center.  Then  1900  volts  are  applied  to  the  wire  from 
a  blaster  to  melt  and  initiate  a  Bpark  at  the  center  of  the 
charge . 

Experimental  Results.  Two  series  of  decoupling  experi¬ 
ments  were  performed  In  plaster  of  parls.  In  each  series  the 
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weight  of  the  charge  was  kept  constant.  In  the  first  series 
a  55  nun  charge  was  used,  and  in  the  second  series  a  100  mm 
charge  was  used.  The  distance  between  the  center  of  the 
charge  and  the  receiver  is  kept  constant.  The  source  detec¬ 
tor  distance  is  15*3  centimeters.  The  radial  components  of 
motion  were  recorded  at  the  surface  of  the  block.  One  re¬ 
ceiver  is  positioned  parallel  to  the  equatorial  plane  of  the 
cavity  and  another  one  perpendicular  to  it. 

For  55  mm  series  the  cavity  radius  is  varied  from  0.47 
centimeter  to  3.2  centimeters,  and  for  100  mm  series  he 
cavity  radius  is  varied  from  0.63  centimeter  to  3.2  .‘.mtime- 
ters.  The  charges  diameters  are  0.47  centimeter  and  0.6  cen¬ 
timeter  . 

Typically,  the  tamped  shots  were  made  in  the  following 
manner.  The  charge  was  placed  at  the  center  of  the  cylindri¬ 
cal  core  during  the  casting.  To  protect  the  charge  from  mois¬ 
ture  the  charge  was  painted  with  oil  paint  before  casting  in¬ 
side  the  specimen. 

The  tamped  shot  gave  small  seismic  signals  in  compari¬ 
son  with  the  cavity  shots.  Post  shot  drilling  of  the  hole 
showed  a  crushing  zone  around  the  charge.  The  diameter  of  the 
cavity  left  after  the  shot  is  as  follows:  for  55  ram  charge, 

7  millimeters;  for  100  mm  charge, 1.1  millimeters.  This  indi¬ 
cates  that  some  of  the  energy  released  by  the  explosions  was 
spent  in  crushing  the  material.  The  cracking  zone  cannot  be 
be  recognized.  Hence,  the  energy  released  by  the  explosion 
was  used  in  an  irreversible  process  leaving  only  part  of  it 
available  for  the  radiated  seismic  energy. 
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Nlcholls  and  Duvall  have  Investigated  both  theoretically 
and  experimentally  the  effect  of  the  ratio  characteristic  im¬ 
pedance  on  explosion  generated  strain  amplitudes  and  strain 
energies.  Although  their  theory  postulated  is  for  plane 
waves,  observed  data  agree  with  predictions  based  on  shock 
wave  theory.  According  to  their  results,  the  stress  in  the 
medium  may  be  calculated  from 


Medium  Stress 


detonation  pressure. 


Where  R  is  the  ratio  of  characteristic  Impedance  of  explosive 
to  rock  and  N  is  determined  to  be  five  from  the  experimental 
data.  N  may  be  defined  as  the  ratio  of  mass  flowing  through 
the  reflected  shock  front  per  unit  area  in  unit  time  to  the 
mass  flowing'  through  the  incident  shock  front  per  unit  area 
in  unit  time  or  as  the  ratio  of  the  effective  characteristic 
impedance  of  the  reflected  shock  wave  to  that  of  incident 
shock  wave.  If  these  ratios  are  equal,  N~l,  and  the  above 
equation  reduces  to  the  acoustical  law  for  the  transmission 
of  a  wave  across  a  boundary  at  normal  incidence. 

The  detonation  pressure  of  the  explosive  used  in 
the  experiment  (silver  acetylide)  was  calculated  from  the 
following  equation  (Ref.  l). 


p  A  0-1*?  . 
l  +  0.8A 


(l) 


where  P  ■  detonation  pressure  in  dyne  per  square  centimeter, 

average  loading  density  in  gram  per  cubic  centimeter, 
D  »  detonation  velocity  in  centimeter  per  second. 
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For  the  silver  acetylide  employed  in  the  experiment,  the  load- 
density  is  0.9  gm/cm^  and  detonation  velocity  1.3  km/sec.  It 
has  to  be  remembered  that  the  detonation  velocity  is  a  func¬ 
tion  of  a  loading  density.  Using  the  above  equation,  the  de¬ 
tonation  pressure  is  given  by, 

P  =  4.0  kilobar. 

The  plaster  of  paris  has  a  density  of  1.1  gram  per  c'j:Dic 
centimeter;  a  longitudinal  propagation  velocity,  c,  of  3.37 
km/sec.,  and  a  characteristic  impedance  (  po)  of  3.63  X  10^ 
dyne  sec/cm^.  The  characteristic  impedance  (  ^  D)  of  the  ex¬ 
plosive  used  is  1  17  X  105  dyne  sec/cm^.  The  impedance  ratio 
is  simply  the  characteristic  impedance  of  the  explosive,  di¬ 
vided  by  the , characteristic  impedance  of  rock,  and  is  given  by 
R  *  0.322.  Using  equation  (l),  the  stress  transmitted  to  the 
medium  is  found  to  equal  to  9.2  kilobar. 

The  compressive  and  tensile  strengths  of  the  plaster  mea¬ 
sured  in  static  failure  are  4l4  and  40  bar  respectively.  The 
dynamic  compressive  and  tensile  strengths  of  the  material  are 
about  twice  the  static  strength. 

Comparing  the  stress  transmitted  to  the  medium  from  a 
tamped  shot  v.'lth  the  dynamic  strengths  of  the  material,  the 
input  stress  to  the  medium  exceeded  the  strength  of  the  mater¬ 
ial.  The  input  stress  is  about  one  order  of  magnitude  larger 
vhan  compressive  strength  and  two  orders  of  magnitude  larger 
than  the  tensile  strength.  Hence,  the  crushed  and  cracked 
zones  should  exist  around  the  tamped  shot. 
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The  variation  of  peak  amplitude  with  cavity  diameter 
from  two  series  of  decoupling  experiments  is  shown  in  Pig.  3 
and  Pig.  4.  Comparing  the  amplitude  of  the  signal  obtained 
from  the  smallest  cavity  (radius  »  0.47  cm  for  55  mgm  charge 
and  0.63  cm  for  100  mgm  charge)  with  the  signal  obtained 
from  the  largest  cavity  (radius  =  3.2  cm),  the  amplitude  of 
the  smallest  cavity  is  about  twice  the  largest  cavity.  This 
is  the  same  for  both  charge  sizes.  Therefore  the  reduction 
in  amplitude  is  independent  of  the  charge  size. 

By  plotting  the  peak  amplitude  as  a  function  of  the 
radius  of  the  cavity  on  logarithmic  paper,  a  straight  line 
can  be  drawn  that  fits  the  observed  points.  The  least  square 
method  is  used  to  get  the  slopes  of  the  lines.  The  slopes  of 
the  lines  obtained  for  55  mgm  charge  and  100  mgm  charge  are 
-0.46  and  -0.55,  respectively.  The  slopes  indicate  that  the 
peak  amplitude  is  inversely  proportional  to  the  cavity  radius 
raised  to  a  power  less  than  one. 

The  maximum  value  of  amplitudes  in  Pigs.  3  and  4  corre¬ 
spond  to  cavity  radii  0.47  cm  and  O.63  cm.  In  the  regions  to 
the  left  of  these  radii  the  amplitudes  start  decreasing  with 
decreasing  cavity  radius.  In  this  region  the  pressure  pro¬ 
duced  by  the  explosion  exceeds  both  the  compressive  and  tensile 
strength  of  the  material.  For  cavity  radii  larger  than  those 
given  above  the  amplitude  decreased  as  the  radius  increased. 

The  pressure  on  the  surface  cf  the  cavity  is  within  the  com¬ 
pressive  strength  of  the  plaster,  but  exceeded  the  tensile 
strength  of  the  plaster.  Therefore  a  cracked  zone  exists  in 
the  medium  surrounding  the  cavity. 
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To  get  a  quantitative  value  of  the  pressure  applied  on 
the  surface  of  the  cavity,  a  shock  wave  theory  would  be  used 
to  calculate  the  pressure  decay  Inside  the  cavity.  To  com¬ 
pare  this  wit*  the  pressure  transmitted  to  the  medium,  the 
amplitude  and  phase  of  the  observed  signal  at  the  surface  of 
the  block  would  be  corrected  for  the  instrument,  propagation 
path  and  geometrical  divergence.  Then  from  the  reduced  sig¬ 
nal  the  transmitted  pressure  would  be  calculated. 
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7.  Determination  of  Depth  of  an  Earthquake  Focus 
from  the  Spectrum  of  Body  Waves. 

Statement  of  the  Problem.  It  is  known  that  the  spectrum 
of  the  seismic  body  waves  as  recorded  at  a  seismographic  sta¬ 
tion  is  influenced  by  several  factors  which  include  the  dept  a 
of  the  earthquake  focus.  The  purpose  of  this  investigation 
is  to  find  a  method  wnich  will  enable  one  to  infer  and  hope¬ 
fully  determine  approximate  focal  depth  from  the  spectrum  of 
seismic  body  waves  as  recorded  at  a  single  or  a  group  of  seis* 
mographlc  stations. 

Theoretical  Background  of  the  Work.  Much  of  the  work  cf 
this  investigation  is  based  on  the  paper  by  Haskell  (1953)  in 
which  he  developed  a  matrix  relation  between  the  motion  stress 
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vector  (I.e.  the  horizontal  and  vertical  components  of  par¬ 
ticle  velocity  and  the  normal  and  tangential  c  ..jponents  of 
stress)  at  an  interface  of  multilayered  media  and  the  ampli¬ 
tudes  of  plane  waves  below  this  interface.  By  an  extension 
of  Haskell's  work,  a  matrix  relation  was  formulated  which 
enables  one  to  calculate  the  amplitudes  of  plane  waves  in 
the  layers  of  the  source  crust  from  the  spectrum  of  body  waves 
recorded  at  a  selsmographic  station. 

Sirce  the  general  matrix  relations  developed  by  Haskell 
are  familiar,  the  theory  will  not  be  presented  in  this  sum¬ 
mary.  Suffice  it  to  note  that  from  the  record  of  dilatational 
waves  at  a  selsmographic  station  one  can  calculate  the  fre¬ 
quency  spectrum  of  Uq/c  and  Wq/c.  Using  the  matrix  relations, 
one  can  obtain  Anv  ,  the  wave  potential  of  ar.  incident  P 
wave  at  the  receiver  crust.  By  an  extension  of  the  matrix  de¬ 
velopment  one  can  then  obtain  the  amplitudes  of  the  plane 
waves,  i.e.  A' (p)»  A'(pb  U(P)*  and  uo(P)»  in  the  sourcts 
layers  which  have  given  rise  to  dilatational  waves  in  the  man¬ 
tle. 


To  obtain  focal  depth  from  these  derived  wave  potentials 

two  methods  were  considered. 

First  method.  If  one  calculates  both  the  dilatational 
and  rotational  wave  potential  distribution  of  particular 
frequency  at  all  points  situated  in  a  vertical  column 
through  the  epicenter  from  the  P  waves  recorded  at  a 
seipmographic  station  i  and  sums  the  derived  wave  poten¬ 
tial  distribution  with  respect  to  i  (i.e.  one  obtains  the 
Bum  of  the  wave  potential  distribution  derived  from  all 
selsmographic  stations  around  the  epicenter),  one  can 
obtain  the  potential  distribution  in  the  neighborhood 
of  the  focus.  If  a  source  is  present  in  the  layered 
crust,  theoretically  one  should  obtain  a  maxima  in  the 
vicinity  of  the  source.  Therefore,  the  point  at  which 
the  computed  potential  distribution  shot*  a  Laxiroa  will 


78 


roughly  correspond  to  the  position  of  the  focus  In 
the  crust. 

There  are  two  major  disadvantages  which  render 
the  application  of  this  method  difficult.  They  ares 

(a)  Record'1  from  a  large  number  of  seismographic 
stations  need  to  be  analyzed  and 

(b)  The  travel  time  from  the  focus  to  the  seismo- 
graphic  station  needs  to  be  known  very  accurately. 

The  wave  potential  distribution  should  dhow  a 
maxima  near  the  focal  region  provided  that  the 
potential  distribution  at  one  particular  instant 
is  considered.  In  order  to  obtain  the  distribu¬ 
tion  of  the  wave  potential  at  one  particular  in¬ 
stant  in  the  source  crust,  one  needs  to  consider 
both  the  amplitudes  and  phases  of  the  wave  poten¬ 
tials  A’(P)i  A!'(p)>  oJ(P)>  and  wi'(P).  Since 
the  time  of  transmission  from  the  focus  to  the 
recording  station  affects  the  phases  of  these  wave 
potentials,  the  transmission  time  needs  to  be 
known  very  accurately  for  all  the  recording  sta¬ 
tions.  (Although,  inaccuracy  in  the  determina¬ 
tion  of  the  origin  time  of  the  earthquake  will 
not  matter,  since  this  will  affect  the  phases  of 
the  wave  potentials  for  all  the  stations  equally.) 

Because  of  these  difficulties,  this  method  has  not  been 
tried  yet. 

Second  method.  The  waveform  recorded  at  a  seismographic 
station  is  a  transient  time  function.  By  a  matrix  formu¬ 
lation  described  earlier,  one  is  able  to  calculate  the 
wave  potentials  in  the  layers  of  the  source  crust  from 
the  recorded  waveform.  These  wave  potentials  include 
the  wave  potentials  of  the  upgoing  and  downgoing  P  waves 
and  those  of  the  upgoing  and  downgoing  S  waves  in  all 
the  layers  of  the  source  crust.  At  any  point  in  the 
source  crust  the  waveform  due  to  any  of  these  potentials 
is  a  transient  time  function.  However,  the  time  of  ini¬ 
tiation  of  these  four  pulses  should  not  be  the  same;  it 
is  only  at  the  focus  that  all  the  four  pulses  should 
start  at  the  same  Instant .  If  one  calculates  the  wave¬ 
forms  due  to  these  four  waves  at  discrete  points  in  a 
vertical  column  through  the  epicenter ,  theoretically 
only  at  one  point  jhould  the  four  transients  initiate 
at  the  same  instant,  and  that  point  should  correspond 
to  the  focus. 

To  test  how  this  method  works,  an  earthquake  of 
magnitude  7  (depth  of  focus  25  Km)  with  epicenter  near 
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Haiti  was  chosen  and  seismograms  recorded  at  Golden, 
Colorado,  were  analyzed  using  a  standard  crustal  model. 
The  spectrum  of  the  derived  wave  potentials  indicated 
a  focus  at  the  crust-mantle  boundary.  It  was  at  that 
state  the  underlying  principle  of  the  second  method  was 
re-examined.  It  appeared  that  the  method  would  not 
work.  At  the  present  time,  we  are  inclined  to  think 
that  the  chances  that  one  will  be  able  to  obtain  focal 
depth  uniquely  from  observations  at  one  single  station 
are  rather  dim,  since,-  as  we  shall  see  later,  tie  spec¬ 
trum  of  the  body  waves  is  influenced  by  factors  other 
than  the  depth  of  earthquake  focus. 

Present  Development.  We  also  considered  the  question  why 
the  spectrum  of  the  seismic  body  waves  be  at  all  Influenced 
by  the  depth  of  the  earthquake  focus  and  the  conclusion  we 
reached  is  that  the  following  factors,  besides  the  character¬ 
istics  of  the  recording  instruments,  will  influence  the  spec¬ 
trum: 

(a)  source  spectrum  and  the  dimension  of  the  source  - 
this  includes  (l)  the  (time)  wave  shape  of  the 
stress  field  at  various  parts  of  an  extended  source 
and  (2)  the  effect  of  integration  of  waveforms 
radiated  from  different  parts  of  the  source  as  the 
recording  station 

b)  radiation  pattern  at  the  source  region 
the  leering  of  the  source  crust 
’d)  the  position  of  the  source  in  the  layered  crust. 

To  demonstrate  the  effect  of  factors  (b),  (c)  and  (d), 

let  us  consider  a  point  source  at  A  in  the  homoger.eous  crust 

shown  below. 
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If  the  reflection  coefficient  at  the  free  surface  is  -1, 
pP  will  be  180  degrees  out  of  phase  with  respect  to  P  and  it 
will  be  delayed  by  a  time  T.  At  the  receiving  station  one 
would  observe  a  waveform  which  is  the  sum  of  pP  and  P. 

Let  P  wave  denoted  as  a  time  function  f^(t)  having  fre¬ 
quency  spectrum  F^iu)  and  pP  be  denoted  as  fg(t)  having  fre¬ 
quency  spectrum  Fg(«ju).  Then  the  frequency  spectrum  of  the 
sum  of  the  two  waveforms  will  be 


p(u) 
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Assuming  that  the  radiation  along  the  directions  AB  and 
AE  are  the  same,  we  have  f2(t)  the  same  as  f^t),  except  that 
f2(t)  is  180  degrees  out  of  phase  and  delayed  by  a  time  T  with 
respect  to  f^(t).  Hence, 

F2M  -  F1(w)  exp[  -  i(ir+u>T5\ 


Therefore, 

fx(t)  +  f2(t)  -  -L.  e1"6  PjW  [l  +  «'  1(T  +  loT>]  <!<*. 
Hence, 

PM  =  F^w)  ^1  +  e’  i(1i  +WT)] 

Therefore,  the  effect  of  summing  the  two  waveforms  is  the  same 
as  filtering  f^(t)  by  a  filter  whose  frequency  response  func¬ 
tion  is  given  by 


KM  -  1  +  e"  1(tr  +WT) 
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Amplitude  response  function  of  this  filter  is  of  the  form 
shown  below 


tt  iTr  tv  +tt 
toT” 


Thus  one  can  see  that  certain  discrete  frequencies  are 
totally  eliminated  in  the  resultant  waveform.  However,  the 
picture  as  shown  above  is  very  simplified.  In  reality  the 
filtering  effect  will  not  be  as  shown  in  the  figure  because 
of  the  following  factors: 

(1)  because  of  influence  of  direction  on  radiation 
pattern,  the  strength  of  the  Bignal  propagated 
along  the  directions  AB  and  AE  are  not  the  same 

(2)  the  laye  ring  of  the  source  crust  will  influence  the 
relative  strengths  of  pP  and  P;  the  layering  will 
introduce  internal  reflections  which  will  also  be 
observed  on  the  seismogram 

(3)  the  reflection  coefficient  at  the  free  surface, 
in  general,  1b  not  -1. 

Because  of  the  above  factors,  none  of  the  frequency  com¬ 
ponents  will  be  completely  eliminated;  however,  there  will  be 
minima  in  the  frequency  spectrum  at  certain  discrete  frequen¬ 
cies.  These  characteristics  of  the  frequency  spectrum  could 
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be  utilized  to  determine  focal  depth.  The  procedure  which 
has  been  recently  thought  about  is  outlined  below: 

(1)  assume  (i)  an  impulsive  point  source  at  some  depth 
in  the  layered  crust,  (ii)  radiation  is  same  in  all 
directions,  and  (iii)  some  reasonable  value  of  re¬ 
flection  coefficient  at  the  free  surface 

(2)  calculate  the  dilatational  waves  emerging  at  the 
base  of  the  crust  at  an  apparent  velocity  c.  This 
time  function  is  the  weighting  function  of  the 
filter  which  results  from  compositing  the  P  waves 
with  pP,  sP  and  internal  reflections  in  the  source 
layers 

(3)  by  process  of  inverse  filtering,  remove  the  effect 
of  filtering  described  earlier  and  *hus  obtain  the 
frequency  spectrum  of  the  source  which  would  have 

'  been  observed  at  the  focus  had  there  been  no  reflec¬ 

tors,  such  as  free  surface  and  internal  interfaces. 

If  the  assumed  focal  depth,  radiation  pattern,  reflection 
coefficient  at  the  free  surface  are  correctly  chosen,  one 
should  obtain  a  reasonably  smooth  frequency  spectrum.  On  the 
other  hand,  if  these  parameters,  particularly  the  focal  depth, 
are  chosen  erroneouply,  one  ehould  obtain  maxima  and  minima 
in  the  spectrum  at  certain  discrete  frequencies.  Thus  by  assum¬ 
ing  focus  at  different  depths  in  the  layered  crust,  by  a  method 
of  trial  and  error,  one  shou'd  be  able  to  locate  the  roost  prob¬ 
able  position  of  the  focus  by  the  criteria  of  the  smoothness 
of  the  calculated  frequency  spectrum. 

It  should  be  pointed  out  that  the  proposed  method  cannot 
uniquely  determine  the  focal  depth.  Also,  in  the  light  of  the 
foregoing  discussions,  it  appears  that  the  chances  that  one 
will  be  able  to  find  a  method  which  can  uniquely  pinpoint  the 
earthquake  focus  from  observations  at  one  single  station  are 
rather  dim.  The  reason  for  this  lack  of  uniqueness  is  that 
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the  spectrum  of  the  body  waves  is  influenced  by  several  fac¬ 
tors  other  than  the  focal  depth  (examples  of  these  factors: 
radiation  pattern,  spectrum  at  the  source  region,  etc.)  and 
from  observations  at  one  station  it  is  not  possible  to  elim¬ 
inate  the  effect  of  other  factors.  However,  the  method  outlined 
above  could  give  the  most  probable  position  of  the  earthquake 
focus . 


Reference 

Haskell,  N.A.  (1953).  The  dispersion  of  surface  waves  on 
multilayered  media.  Bull.  Seism,  Soc.  Am.  43,  17-34. 

8.  Factors  Affecting  the  Determination  of 
Phase  Velocities. 

Phase  velocity  curves  obtained  from  seismograms  are  not 
always  smooth,  nor  do  they  have  the  expected  shape  of  theoreti¬ 
cal  phase  velocity  curves.  As  a  result  cf  difficulties  en¬ 
countered  in  the  determination  and  interpretation  of  these 
curves  considerable  time  has  been  spent  during  the  last  two 
quarters  of  this  year  in  examining  into  the  causes  of  these 
irregularities.  The  investigations  are  here  summarized. 

In  particular,  since  the  method  of  processing  of  data 
may  contribute  notably  to  the  irregularities,  we  have  examined 
quantitatively  the  effects  of  factors  such  as  ramp  functions 
and  data  windows  on  the  phase  velocity  curves.  For  this  pur¬ 
pose  nine  pairs  of  stations  were  selected  more  or  less  at  ran¬ 
dom  from  sets  of  Rayleigh  wave  data  for  four  earthquakes.  The 
identical  digitized  data  for  each  station  pair  were  processed 
ty  different  methods  and  the  results  compared.  The  earthquakes 
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and  station  pairs  examined  pre  given  in  Table  1. 

In  normal  procedures  records  are  digitized  at  approximate¬ 
ly  one  second  intervals.  The  digitized  data  a-e  then  filtered, 
smoothed  and  decimated  using  a  six  point  filter  (see  Hamming, 
1962,  p  317)  and  the  trend  and  mean  removed.  A  Fourier  analy¬ 
sis  is  then  performed  by  applying  Simpson1 s  integration  tech¬ 
nique.  Corrections  for  instrument  magnification  and  phase  re¬ 
sponse  are  made,  and,  finally,  the  phase  velocity  computations 
are  performed  by  Sato’s  method  ( 1955 ) • 

Except  for  a  deconvolution  procedure,  variations  in  the 
data  processing  methods  (  that  is,  the  way  in  which  the  digi¬ 
tized  data  are  filtered  and  smoothed,  the  way  in  which  the 
slope  of  the  ramp  function  is  determined,  and  the  stage  at 
which  the  correction  for  the  mean  is  to  be  made)  take  place 
before  the  Fourier  analysis. 

Seven  methods  of  data  processing  preceding  phase  velocity 
computations  were  applied. 

1.  The  ramp  angle  was  determined  by  means  of  a  least 
square  fit  of  the  set  of  digitized  data  to  a 
straight  line.  By  correcting  for  the  slope  and 
intercept  of  the  straight  line  the  trend  and  the 
mean  are  removed  from  the  data.  Processing  of  the 
data  then  proceeded  as  indicated  above. 

2.  The  ramp  angle  was  determined  by  computing  the 
average  of  the  first  half  of  the  digitized  data 
series  and  the  average  of  the  second  half,  then 
passing  a  straight  line  thru  the  two  points  so 
located  at  one-fourth  and  three-fourths  of  the 
time  duration  of  the  data.  The  filtered  data 
were  corrected  for  this  ramp  angle,  and  then  the 
correction  for  the  mean  made  by  subtracting  the 
difference  between  the  mean  of  the  filtered  and 
ramp-corrected  data. 
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TABLE  I 


08  August  1963,  Pox  Is.  EBMT 

FLO 

21  August  1963,  Caribbean  Is.  ATL  - 

07  September  1963,  Ascension  Is.  ATL  - 

MDS  - 


OIMA 

FRMA 

AAM 

TDNM 


-  BLA 

-  ATL 

FLO 

RTNM 

FRMA 

-  MDS 

-  BLWV 

-  BRPA 

-  DU  OK 


06  November  1963,  New  Guinea 
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3.  The  ramp  angle  was  computed  as  In  process  2,  then 
the  filtered  data  were  corrected  for  the  slope  and 
Intercept  of  the  ramp. 

4.  The  data  were  filtered  and  the  mean  only  was  re¬ 
moved;  no  correction  was  made  for  the  ramp  angle. 

3.  The  ramp  was  determined  as  In  process  2,  then  a 
boxcar  function  with  height  A-(L  -  R)/2  was  Intro¬ 
duced.  Here  A  Is  the  mean  of  the  filtered  data, 

L  Is  the  mean  of  the  first  half  of  the  digitized 
data,  R  the  mean  of  the  second  half.  This  method 
of  computing  the  boxcar  height  was  chosen  so  that 
It  would  also  be  a  function  of  the  amplitudes  of 
the  data. 

6.  The  ramp  was  determined  as  In  process  2.  Data 
processing  continued  as  in  that  process  but  with 
the  application  of  a  hanning  window  to  the  fil¬ 
tered,  and  ramp  and  mean  corrected  data.  The 
hanning  window  Is  given  by 

a(^)  -  0.5  [l  +  coa  f(-— y-Tj]  <  Itjczr  .  1  -  1,2,. ...n 

-  0  ,  |n|>2T 

where  n  is  the  number  of  data  points,  T  is  the 
time  of  the  ith  data  point  and  T  is  one*half  the 
time  duration  of  the  data  (see  Blackman  and  Tukey). 

Y.  A  deconvolution  procedure  was  applied.  This  was 
based  on  the  assumption  that  the  finite  length  sur¬ 
face  wave  train  was  the  product  of  an  infinite  wave 
train  and  a  boxcar  function  in  time.  The  resulting 
spectrum  is  the  convolution  of  the  spectrum  of  the 
infinite  wave  train  with  the  spectrum,  of  the  boxcar 
function.  Deconvolution  to  find  the  spectrum  of 
the  infinite  wave  train  involves  the  solution  of 
integral  equations.  Since  the  convolution  spectrum 
is  that  obtained  by  Fourier  analysis  of  the  data, 
and  since  knowing  the  length  of  the  data  we  can  com¬ 
pute  the  spectrum  of  the  boxcar  function  of  unit 
height,  it  is  possible  to  solve  the  Integral  equa¬ 
tion  for  the  spectrum  of  the  infinite  wave  train. 

In  the  deconvolution  procedure  the  phase  velocity  curves 
obtained  are  in  general  jagged,  with  deep  minima.  There  is 
little,  if  any,  resemblance  to  the  shape  of  theoretical  phase 
velocity  curves.  From  this  it  appears,  that  the  assumption  of 


multiplication  by  a  boxcar  function  in  time  may  be  invalid. 
This  implies  that  the  physical  situation  of  the  data  them¬ 
selves  corresponds  to  a  signal  of  finite  time  duration. 

The  application  of  a  harming  window  also  gave  erratic 
phase  velocities.  The  Fourier  amplitude  spectra  have  the 
same  shape  as  without  hanning,  but  are  reduced  in  amplitude 
and  are  smoother.  The  Fourier  phase  spectra  are  generally 
quite  different  from  the  phase  spectra  found  witho\t  harming. 

In  short,  neither  deconvolution  nor  the  application  of  a 
, harming  window  gave  favorable  results.  No  further  discussion 
will  be  given  of  these  two  procedures. 

There  remain  procedures  l  through  5.  To  examine  the 
various  parameters  of  the  data  processing  procedures  for 
their  effect,  on  phase  velocity  determination,  computations 
from  four  representative  pairs  of  stations  from  those  given 
in  Table  1  are  here  presented. 

1.  Pair  ATL-ELO,  21  August  1963  (see  Fie.  1*) 


Method  1.  For  ATL  the  least  squares  ramp  angle,©, 
ii  +0.1°  and  h/V*  ■  0,001,  where  h  is  the  value 
of  the  mean  and  (I- is  the  standard  deviation. 

For  FLO,  ©  «  0.0  and  C 'T •  0.000.  Obviously  this 
is  a  very  good  set  of  data. 

Method  2.  For  “"thod  ©  ■  +2.2°  and  0.2°,  respec¬ 
tively,  whi*  ..  »  o.  This  gives  exactly  the 
same  phase  velocities  as  method  1. 


In  fact  methods  2  through  5  give  good  agreement 
with  the  least  squares  method.  The  maximum  ramp 


*ln  the  figures  the  entire  phase velocity  curves  from  8  to 
72  seconds  are  presented  for  purposes  of  comparison  of  the 
various  processing  methods.  It  should  be  noted  th?it  in  some 
cases  the  curves  should  be  terminated  at  periods  correspond¬ 
ing  to  minima  of  the  Fourier  amplitude  spectra.  The  ampli¬ 
tude  spectra  are  indicated  in  the  figures. 
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angle  Is  2.2°  and  the  maximum  value  of  h /tf  is 

0.081. 

The  probable  reason  for  this  is  the  relative¬ 
ly  low  amplitudes  of  the  surface  waves.  The 
possible  range  of  the  ramp  angle  and  the  value 
of  the  mean  are  much  restrained. 

2.  Pair  ATL-RTNM,  07  September  1963  (see  Pig.  2). 

Method  1.  The  curve  Is  net  so  smogth  as  fur  the  pre¬ 
ceding  pair.  For  ATL,  9  «  -2.5  and  h/cr*  «  0.003; 
for  RTNM  9  «  0.0  and  h/p-  »  .25. 

Method  2.  This  method  gives  almost  the  same  phase 
velocity  curve  as  that  from  method  1,  but  smooth¬ 
er.  For  ATL  9  *  -5.8°;  for  RTNM,  9  -  +4.1°.  The 
mean  is  zero.  Evidently  for  the  period  range 
8-72  seconds  the  effects  of  the  ramp  angles  can¬ 
cel  and  since  the  mean  is  zero  a  smoother  phase 
velocity  is  obtained. 

Methods  3  and  4.  These  methods  each  have  h/cr  »  0.0 
and  yield  the  same  angles  as  above,  thus  giving 
the  same  phase  velocities. 

Method  5.  Here  tyAr  *  0.01  for  ATL  and  0.236  for 
RTNM.  The  resulting  phase  velocity  curve  oscil¬ 
lates  about  the  curve  found  by  method  1.  This 
may  be  caused  by  the  different  ramp  angles. 

3.  Pair  FRMA-BLWV,  06  November  1963  (see  Fig.  3). 

Method  1.  For  FRMA,  9  -  -0.1°,  and  for  BLWV  9  »  0.1°, 

h/(T=  0.012,  and  for  BLWV  9  -  +0.1°,  h/£r*  =  0.022. 

Method  2.  The  phase  velocities  match  those  of  method 
1.  But  for  FRMA  9  =  3.7°  and  for  BLWV  9  «  22.2°, 
with  zero  means.  The  period  range  of  interest 
(8-72  sec.)  and  the  relative  time  durations  of  the 
data  (FRMA,  579  sec;  BLWV,  327  sec.)  must  account 
for  the  matching  of  the  phase  velocities.  The 
angle  of  22.2°  would  seem  to  be  too  large  to  allow 
reasonable  results. 

Method  3.  The  values  of  9  are  -0.8°  and  -3.9°#  re¬ 
spectively;  the  means  are  zero.  The  phase  veloci¬ 
ties  are  very  nearly  identical  to  those  of  method  1. 

Method  4.  Pnase  velocity  curve  is  irregular. 

Method  5.  9  ■  3-7°  and  22.2°,  respectively,  and 

h/tf"  ■  0,138  and  0,064.  The  respective  means 
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are  -29.9  and  +30.2.  It  appears  that  the  effect 
of  boxcar  functions  approximately  equal  in  magni¬ 
tude  but  opposite  in  sign  cancel  for  the  period 
range  8-72  sec. 

4.  Pair  EBMT-BLA,  08  August  1963  (see  Fig.  4). 

Method  1.  For  EBMT  6  -  -7.2°  and  h/f  -  0.014;  for 
BLA  0  ■  -0.4°  and  h/j-  =  0.34,  The  combination 
of  the  large  6  at  the  first  station  and  large 
h/cf  at  the  second  probably  account  for  the  appear¬ 
ance  of  the  curve. 

Methods  2-4.  These  methods  give  approximately  the 
same  type  of  oscillating  curves.  For  all  three 
methods  ramp  angles  are  -21.2°  and  27  0°  with 
zero  means. 

Method  5.  Ramp  angles  the  same  as  in  2-4,  but  h Jqr 
*  0.027  and  0.182,  respectively. 

Although  the  data  from  these  pairs  of  stations  are  few, 
it  Is  possible  to  propose  some  tentative  results.  If  the 
maximum  period  of  interest  is  approximately  one-third  or  less 
of  the  time  duration  of  the  data,  the  effects  of  boxcar  func¬ 
tions  of  about  the  same  magnitude  but  of  opposite  sign  at  each 
of  the  two  stations  cancel  when  computing  phase  velocities. 

The  same  seemc  to  be  true  of  the  ramp  angle.  It  also  appears 
that  a  boxcar  of  seme  relatively  small  amplitude  and/or  a 
small  ramp  angle  may  be  allowable  but  the  relation  of  the 
magnitudes  of  the  two  functions  to  other  parameters  is  ambigu¬ 
ous  . 


Relatively  low  araplil  'e  surface  waves  give  the  better 
results  for  the  first  five  data  processing  methods.  This  is 
attributable  to  the  restriction  in  the  range  of  the  mean  and 
of  the  ramp  angle. 

The  procedure  of  using  the  first  half  and  the  second  half 
of  the  data  to  compute  the  ramp  angle  (to  be  called  the  half 
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and  half  method)  always  gives  angles  larger  than  the  least 
squares  method.  Filtering  of  the  digitized  data  usually  in¬ 
troduces  a  ramp  angle  large  enough  that  Fubtracting  the  slope 
of  the  original  data  from  the  filtered  data  there  still  re¬ 
mains  a  ramp  angle  that  is  two  to  four  times  the  ramp  angle 
in  the  original  (digitized)  data.  Thus  the  data  that  are 
Fourier  analyzed  often  have  a  ramp  angle  larger  than  that  in 
the  original  data.  Note  that  these  ramp  angles  are  computed 
by  the  half  and  half  method. 

The  ramp  argle  of  the  original  data  computed  by  the 
least  squares  method  is  always  less  than  the  ramp  angle  com¬ 
puted  by  the  half  and  half  i^thod.  With  the  least  squares 
method,  the  ramp  angles  of  the  original  and  filtered  data  are 

within  0.2°  .of  each  other,  so  that  the  ramp  angle  of  the  data 

o 

which  are  Fourier  analyzed  is  less  than  +1,0  in  the  majority 

o 

of  cases,  and  often  less  than  +  0.4  . 

While  the  least  squares  method  of  computing  the  ramp 
angle  was  taken  as  the  standard  method  against  which  the 
others  were  compared,  this  does  not  necessarily  mean  it  is 
the  best  method.  The  half  and  half  method  often  gave  phase 
velocity  curves  whose  shapes  were  very  nearly  the  s**  je  as 
the  curves  from  the  least  squares  method,  even  though  they 
did  not  agree  in  numerical  value.  Part  of  the  reason  for 
this  is  apparently  compensating  effects  in  each  of  the  two 
eets  of  data.  It  may  be  that  the  half  and  half  method  gives 
a  better  representation  of  the  base  line  thru  the  data  be¬ 
cause  of  its  dependence  upon  the  means  of  the  first  half  and 
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and  second  half  of  the  data  aeries,  whereas  the  least  squares 
method  is  related  to  the  mean  of  the  entire  data  series .  This 
comparison  is  complicated  by  the  effect  that  the  half  and  half 
method  leaves  large  ramp  angles  in  the  data. 

Another,  and  possibly  more  serious,  problem  is  that  of 
over  correction.  For  some  data  sets  that  have  small  ramp 
angle,  it  was  attempted  to  again  remove  the  mean  and  ramp 
angle  by  the  least  squares  method  before  Fourier  analyzing. 
But,  unexpectedly,  the  ramp  angle  increases,  though  the  mean 
becomes  closer  to  zero. 

Intuitively,  the  least  squares  method  would  seem  to  be 
the  bp*t  method  for  removing  the  ramp  function  from  the  data. 
Also,  the  ramp  angles  computed  by  the  least  squares  technique 
are  generally  smaller,  and  therefore  more  acceptable,  than 
the  angles  computed  by  the  half  and  half  technique.  However, 
since  neither  of  the  methods  is  able  to  completely  remove  the 
ramp  angle,  more  work  in  this  area  is  required. 
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Work  already  published  or  submitted  for  publics ticn. 

Work  mentioned  in  this  section  has  already  been  reported 

in  detail.  Consequently,  by  way  of  summary,  abstracts  of  the 

papers  in  question  are  here  presented. 

1.  "The  Focal  Mechanism  of  the  Alaska  Earthquake  of  28  March 
1964  and  of  Its  Aftershock  Sequence,"  William  Stauder 
and  G.  A.  Bollinger.  Jour.  Geoph-a.  Res.  71*  (1966): 
5283-5296. 


Abstract 

Focal  Mechanisms  have  been  determined  for  one  pre¬ 
shock,  for  the  main  shock,  and  for  more  thc.i  twenty-five 
aftershocks  of  the  Alaska  earthquake  of  28  March  1964. 

For  the  main  shock  a  single  nodal  plane  with  strike  N66°E, 
dip  850  southeast  is  determinable  from  the  polarity  of 
the  P  wave.  This  plane  may  be  taken  either  as  a  plane 
normal  to  the  fault  motion  (thrust  faulting)  or  as  the 
fault  plane  (dip-slip  motion  on  a  near  vertical  plane). 

A  combination  of  P  wave  first  motion  and  S  wave  polariza¬ 
tion  data  make  possible  the  determination  of  both  nodal 
planes  in  each  shock  studied  of  the  aftershock  sequence. 
One  of  these  planes  is  near  vertical  and  resembles  cloge- 
ly  the  nodal  plane  of  the  main  shock,  the  other  dips  5 
to  15°  to  the  northwest  or  north.  For  earthquakes  of  the 
Kodiak  Island  region,  the  near  vertical  plane  has  the 
same  orientation  as  that  of  the  main  shock.  For  earth¬ 
quakes  to  the  east  of  Prince  William  Sound,  this  plane 
shows  a  systematic  change  in  orientation  corresponding 
to  the  change  in  trend  of  the  tectonic  features.  Three 
earthquakes,  with  foci  at  increasing  depths  along  the 
line  of  greatest  flexure  of  the  tectonic  features  differ 
from  the  others.  The  difference  in  character  of  these 
foci  probably  provides  an  important  clue  to  the  right 
Interpretation  of  the  motion  in  the  main  shock.  While 
the  focal  mechanism  solutions  for  the  shocks  here  studied 
are  subject  to  the  same  ambiguity  of  the  interpretation 
as  that  of  the  main  shock,  criteria  which  favor  a  thrust 
hypothesis  are  advanced  from  the  inter-relation  of  the 
foci.  From  dislocation  theory  it  5s  shown  that  differ¬ 
ential  slip  and/or  a  dipping  thrust  plane  explain  satis¬ 
factorily  the  observed  vertical  displacements  at  the 
surface. 
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2.  "The  S-Wave  Project  for  Focal  Mechanism  Studies,  Earth¬ 
quakes  c.-  1963#"  William  Stauder  and  Q.  A.  Bollinger, 
Bull,  Seism.  Soc.  Am.  56  (1966):  1363-1371. 

Abstract 

P  and  S  wave  data  for  thirty-five  earthquakes 
selected  from  among  the  larger  earthquakes  of  1963  have 
been  investigated.  Focal  mechanism  determinations  for 
twenty-six  of  these  shocks  are  here  presented.  The 
solutions  are  based  upon  a  combination  of  a  graphical 
and  a  computer  method  for  determining  the  poles  of  the 
nodal  planes.  In  all  cases  it  has  been  found  that  the 
mechanism  may  be  adequately  represented  by  a  double 
couple  as  an  equivalent  point  source  of  the  focus,  al¬ 
though  in  some  few  instances  a  single  couple  cannot  be 
excluded  as  a  possible  alternate  interpretation.  The 
solution  of  a  mid-Atlantic  earthquake  of  November  17, 

1963  is  presented  as  an  example  of  a  focus  which  clearly 
conforms  to  the  double  couple  model. 

Special  attention  is  called  to  the  solutions  for  a 
series  of  earthquakes  in  the  Kurile  Islands,  and  to  three 
earthquakes  of  the  mid-Atlantic. 


3,  "Master.  Curves  for  the  Response  of  layered  Systems  to 
Compressional  Seismic  Waves."  Luis  M.  _-<:rnandez,  Bull. 
Seism.  Soc,  Am.  56  (1966):  1363-1371. 

Abstract 

The  layers  of  the  earth's  crust  act  as  a  filter 
with  respect  to  seismic  energy  arriving  at  a  given  sta¬ 
tion.  Consequently  the  motion  recorded  at  the  surface 
depends  not  only  on  the  frequency  content  of  the  source 
and  on  the  response  characteristics  of  the  recording 
instrument,  but  also  on  the  elastic  parameters  and  thick 
ne:jses  of  the  transmitting  layers.  This  latter  depend¬ 
ence  is  the  basis  for  a  method  of  investigating  the 
structure  of  the  crust  and  upper  mantle. 

To  facilitate  this  investigation  a  set  of  master 
curves  for  the  transfer  functions  of  the  vertical  and 
horizontal  component  of  longitudinal  waves  and  their 
ratios  is  presented.  The  calculation  of  these  curves 
is  in  terms  of  a  dimensionless  parameter.  This  calcu¬ 
lation  allows  one  to  group  the  curves  corresponding  to 
different  crustal  models  into  families  of  curves.  The 
characteristics  of  these  curves  are  discussed  from  the 
point  of  view  of  their  "periodicity"  in  the  frequency 
domain  and  of  their  amplitude  in  order  to  investigate 
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the  influence  of  the  layer  parameters.  Considerations, 
either  of  constructive  interference  or  of  Fourier  analy¬ 
sis  of  a  pulse  multiply  reflected  within  the  layer  sys¬ 
tem,  reveal  that  the  amplitudes  of  the  transfer  curves 
depend  on  the  velocity  contrasts  at  the  interfaces  of 
the  system.  The  "periodicity"  or  spacing  of  the  peaks 
depends  on  the  time  lags  between  the  first  arrivals  and 
the  arrivals  of  the  different  reverberations.  Closely 
spaced  fluctuations  correspond  to  large  time  lags,  and 
widely  spaced  fluctuations  to  short  time  lags. 


4.  "Model  Investigation  of  Explosions  in  Prestressed  Media," 
W.  H.  Kim  and  Carl  Klsslinger,  Geophysics,  32  (1967):  633- 

651 . 


Abstract 

Seismic  effects  of  explosions  and  rupture  propagation 
in  prestressed  two-dimensional  models  (Plexiglas,  alumi¬ 
num)  as  well  as  anisotropy  produced  by  the  stress  field 
are  investigated.  An  explosion  in  a  prestressed  medium 
releases  a  portion  of  the  stored  strain  energy  by  one  or 
more  of  the  following  mechanisms:  (l)  formation  of  direc¬ 
tional  cracking,  especially  in  brittle  materials,  (2)  re¬ 
lease  of  strain  energy  in  the  elastic  zone  outside  the 
cavity,,  and  (3)  rupture  propagation.  Phenomena  associated 
with  all  of  these  mechanisms  were  observed  in  the  present 
investigation. 

Explosions  in  prestressed  Plexiglas  produce  cracks 
in  preferred  directions,  the  intensity  of  which  Increases 
with  the  applied  stress.  Explosions  in  prestressed  alumi¬ 
num  sheets  do  not  cause  fracturing  but  rather  plastic  de¬ 
formation  about  the  explosion.  Straight  and  branching 
modes  of  moving  cracks  initiated  from  explosions  in  pre¬ 
stressed  Plexiglas  can  be  explained  on  the  basis  of  stress 
distribution  ahead  of  the  crack  tips. 

Observed  radiation  patterns  resulting  from  explo¬ 
sions  in  prestressed  media  indicate  asymmetrical  radia¬ 
tion  fields  which  are  a  direct  consequence  of  strain 
energy  release  for  the  case  of  aluminum  and  by  the  com¬ 
bined  effects  of  directional  cracking  and  energy  release 
in  the  elastic  zone  for  the  case  of  Plexiglas.  Explosions 
in  prestressed  media  generate  shear  waves.  The  observed 
S-wave  magnitude  increases  sharply  with  the  level  of  the 
existing  stress  field  for  a  given  amount  of  strain  energy 
release.  It  is  concluded  that  this  phenomenon  is  attribu¬ 
table  to  the  effective  conversion  of  energy  release  to 
seismic  radiation  at  high  ambient  stress  fields.  In  other 
words,  the  effectiveness  of  S-wave  generation  is  governed 
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by  the  rapidity  with  which  the  existing  strain  energy  is 
released. 

A  definite  anisotropy  effect  was  observed  in  pre¬ 
stressed  models,  but  this  effect  is  not  large  enough  to 
affect  wave  propagation  in  the  range  of  the  tensile  loads 
applied. 


5.  "Ray  Theory  in  an  Anisotropic  Inhomogeneous  Elastic 
Medium,"  N.  J,  Vlaar,  submitted  to  Bull.  Seism.  Soc. 
Am. 


Abstract 

A  ray  theory  is  developed  for  an  inhomogeneous, 
anisotropic  medium,  based  on  the  concept  of  a  wave 
front  as  the  carrier  of  a  discontinuity  in  particle 
velocity.  The  discontinuity  conditions  of  several 
field  quantities  involved  are  formulated  and  serve  to 
cast  the  problem  in  terms  of  a  partial  differential 
equation.  The  characteristics  of  this  equation  pubse- 
quently  are  identified  in  terms  of  rays,  represented 
in  parametric  for”»  space  curves.  The  energy  trans¬ 
port  along  the  ra.*'  .s  formulated  by  means  of  transport 
equations. 

The  theory  is  applied  in  particular  to  the  case 
of  the  transversely  (horizontally)  isotropic,  vertically 
inhomogeneous  medium.  Equations  for  rays  and  travel 
times  are  obtained  in  the  form  of  integrals,  which  are 
suited  for  serving  the  purpose  of  numerical  computa¬ 
tions  . 


6.  "Determination  of  Earthquake  Fault  Parameters  from  Long 
Period  F  Waves,"  0.  A.  Bollinger,  submitted  to  Journ. 
Qeophys.  Res. 


Abstract 

The  theory  for  propagating  earthquake  sources  pre¬ 
dicts  a  "Doppler-type"  effect  in  the  time  duration  of 
the  radiated  body  wave  pulses.  With  the  assumption  of 
a  constant  fracture  velocity  the  time  duration  of  these 
pulses  is  directly  related  to  the  length  of  the  fault 
but  has  different  spatial  patterns  of  variation  for  the 
unilateral  and  bilateral  faulting  cases .  To  relate  this 
time  duration  of  energy  release  to  the  observed  long- 
period  P  waves,  various  pulseforms  are  postulated  whose 
time  history  is  controlled  by  the  fault  paxameters  of 
length  and  fracture  velocity.  These  pulses  are  convolved 
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with  the  impulse  response  of  a  standard  crust-seismograph 
system  to  yield  theoretical  P  wave  seismograms.  A  rela¬ 
tion  is  thereby  established  between  the  time  durations 
of  the  radiated  pulses  and  those  observed  at  the  surface 
of  the  earth.  This  relationship  is  employed  to  correct 
the  times  of  observed  long  period  P  waves  so  as  to  re¬ 
cover  the  original  time-duration  radiation  pattern.  By 
a  minimum  deviation  technique  the  fault  length  and  frac¬ 
ture  velocity  can  be  found  that  best  fits  this  observed 
pattern.  This  procedure  is  shown  to  yield  reasonable 
results  for  four  different  earthquakes  from  three  dif¬ 
ferent  regions. 


7.  "The  Thickness  of  the  Crust  in  Central  U.S.  and  La  Paz, 
Bolivia,  from  the  Spectrum  of  Longitudinal  Seismic  Waves" 
Luis  M.  Fernandez,  submitted  to  Bull.  Seism.  Soc.  Am. 

Abstract 

The  ratio  of  the  vertical  and  horizontal  components 
of  the  P  phase  provides  a  crustal  transfer  function 
which  is  independent  of  the  frequency  content  of  the 
source  and  is  a  function  only  o'*  the  angle  of  incidence 
of  the  rays  and  the  crustal  parameters  of  the  site  where 
the  observation  is  taken. 

This  ratio  of  the  vertical  and  horizontal  spectra 
corresponds  to  the  tangent  of  the  apparent  angle  of  emer¬ 
gence  and  as  such  this  apparent  angle  is  a  function  of 
the  frequency.  Spectra  of  this  type  have  been  obtained 
in  the  Central  U.S.  and  in  La  Paz,  Bolivia,  for  large 
magnitude  teleaeisms .  The  observed  curves  have  been 
compared  with  theoretical  universal  curves  corresponding 
to  one  and  two  layer  models. 

Average  results  for  several  observations  give  a 
crustal  thickness  in  the  Central  U.S.  of  42  kilometers 
and  a  mean  P  velocity  of  6.6  km/sec.  For  the  Bolivian 
Andes  at  La  Paz,  the  crustal  thickness  obtained  is  of 
64  kms  and  the  mean  P  velocity  of  6,7  km/sec.  These 
results  are  in  good  agreement  with  similar  determina¬ 
tions  obtained  by  independent  methods. 


